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ABSTRACT  The  U.S.  Naval  Observatory  is  faced  with  a  growing  problem. 
Heat  rejection  by  large  buildings  constructed  adjacent  to  the  Observatory  is 
distorting  the  properties  of  the  atmosphere.  This  distortion  refracts  the  incoming 
light  in  an  anomalous  manner,  introducing  uncorrectable  errors  into  celestial 
observations. 

A  building  heats  the  atmosphere  through  the  discharge  of  the  heating  and  air 
conditioning  systems  and  also  by  means  of  heat  convected  off  hot  walls  and 
roof.  As  the  air  is  heated,  its  physical  properties  are  changed.  This  heating  is 
heterogeneous,  influenced  by  the  geometry  of  the  buildings  and  by  the  turbu¬ 
lence  of  the  air.  The  perturbed  atmosphere  then  rises  due  to  the  buoyancy  of 
the  heated  air  and  is  swept  downwind  toward  the  Observatory. 

The  energy  transferred  from  the  buildings  to  the  air  and  the  resulting  change 
in  the  properties  of  the  atmosphere  were  examined  theoretically.  Building  heat 
rejection  was  significant,  particularly  in  the  late  afternoon  and  early  evening, 
throughout  the  year.  Heat  rejection  rates  of  5  MW  for  a  large  building  are 
typical.  When  an  entire  street  is  included,  an  anomalous  telescope  error  of  up  to 
0.3  seconds  of  arc  is  predicted. 
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INTRODUCTION 


The  U.S.  Naval  Observatory  is  faced  with  a  growing  problem.  Heat 
rejection  by  large  buildings  constructed  adjacent  to  the  Observatory  is 
distorting  the  physical  properties  of  the  lower  atmosphere.  This  dis¬ 
tortion  refracts  (bends)  the  incoming  light  in  an  anomalous  manner, 
introducing  uncorrectabl e  errors  into  celestial  observations. 

The  distortion  of  the  atmosphere  is  very  small.  However,  a  primary 
mission  of  the  Naval  Observatory  is  to  determine  and  maintain  a  sophis¬ 
ticated  celestial  reference  system  for  such  applications  as  satellite 
navigation,  missile  guidance  systems,  stellar  tracking  systems,  and  sat¬ 
ellite  station  keeping.  Navigational  errors  of  even  fractions  of  a 
second  of  arc  are  not  tolerable.  The  problem  is  discussed  in  detail  in 
Appendix  A. 

A  building  heats  the  atmosphere  in  two  ways.  Energy  is  transferred 
to  the  atmosphere  with  the  discharge  of  heating  and  air  conditioning 
systems  (HVAC).  Buildings  also  heat  the  air  by  convection  heat  transfer 
off  their  walls.  During  the  morning  and  early  afternoon,  the  walls  re¬ 
ceive  more  solar  radiation  than  can  be  conducted  into  the  building  or 
can  be  transferred  to  the  wind,  and  the  walls  heat  up.  In  the  late 
afternoon  and  early  evening,  the  temperature  of  the  air  decreases,  and 
large  temperature  gradients  exist  between  the  still  hot  walls  of  the 
building  and  the  surrounding  atmosphere.  The  stored  energy  is  then  con- 
vected  to  the  air  blowing  over  these  walls  at  a  rate  proportional  to 
this  temperature  gradient  and  to  the  velocity  of  the  wind. 

As  the  air  is  heated,  its  physical  properties  are  changed.  This 
heating  is  heterogeneous,  influenced  by  the  geometry  of  the  buildings 
and  by  the  turbulence  of  the  air.  The  perturbed  atmosphere  rises  due  to 
the  buoyancy  of  the  heated  air  and  is  swept  downwind  (e.g.,  toward  the 
Observatory). 

The  problem  can  be  summarized  as  a  sequence  of  four  elements: 

1.  Buildings  are  heated  by  the  sun  and  internal  energy  consump¬ 
tion,  and,  in  general,  the  amount  of  energy  involved  is 
proportional  to  the  size  of  the  building. 

2.  A  hot  building  emits  "heat,"  and  this,  in  conjunction  with  the 
wind,  changes  the  downwind  structure  of  the  atmosphere. 

3.  A  perturbed  atmosphere  changes  the  direction  of  travel  of 
incoming  (star)  light  in  an  unpredictable  way,  and  this 
changes  observed  (star)  positions. 

4.  An  additional  error  is  thus  introduced  into  the  observations, 
and  this  adversely  affects  the  accuracy  of  sophisticated  navi¬ 
gation  systems. 
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The  solution  to  this  problem  is  to  restrict  the  size  of  buildings 
constructed  upwind  of  the  Observatory.  To  accomplish  this,  the  Navy 
petitioned  the  District  of  Columbia  Zoning  Commission  to  establish  a 
"precinct  district"  in  the  neighborhood  immediately  surrounding  the 
Observatory.  Within  this  district,  the  height  of  future  construction 
would  be  limited. 

The  Zoning  Commission  was  not  to  be  convinced  by  qualitative  argu¬ 
ments,  however.  Numbers  were  required.  They  would  insist  on  knowing: 
How  much  error  is  being  introduced  into  the  observations? 

Property  changes  of  such  small  magnitudes  occurring  over  regions  of 
several  hundred  yards  are  very  difficult  to  measure.  Instead,  the  sever¬ 
ity  of  the  problem  was  ascertained  theoretically.  The  first  two  ele¬ 
ments,  energy  transferred  from  the  buildings  to  the  air  and  the  result¬ 
ing  change  in  the  structure  of  the  atmosphere,  were  examined  by  the 
Naval  Civil  Engineering  Laboratory  (NCEL).  Using  density  as  the  measure 
of  structure,  a  parametric  study  was  conducted,  varying  building  size 
and  location,  number  of  buildings,  wall  construction,  time  of  day,  time 
of  year,  and  the  direction  of  the  wind.  This  report  discusses  the 
methods  used  by  NCEL  to  make  this  study  and  presents  the  results. 


GEOGRAPHY  OF  THE  SITE 

The  Naval  Observatory  is  located  on  a  72-acre  site  in  northwest 
Washington,  DC.  The  grounds  are  circular,  with  a  300-meter  radius,  as 
illustrated  schematically  on  Figure  1. 

The  telescopes  being  threatened,  the  transit  circle  telescope  and 
26-inch  refracting  telescope,  are  located  approximately  in  the  center  of 
the  grounds.  Figures  2  and  3  show  the  inside  and  outside  of  the  transit 
house.  This  is  the  telescope  of  major  concern.  It  is  used  for  the  fun¬ 
damental  star  position  work  conducted  at  the  Observatory.  This  tele¬ 
scope  is  mounted  so  that  it  can  only  observe  along  one  meridian,  i.e., 
the  telescope  can  only  look  to  the  north  and  the  south.  Stars  and 
planets  are  observed  as  they  "transit"  that  meridian.  Weather  permit¬ 
ting,  observations  are  made  24  hours  a  day  throughout  the  year. 

The  26-inch  refracting  telescope  is  the  largest  telescope  at  the 
Naval  Observatory.  It  is  primarily  used  for  observing  multiple  star 
systems.  Figure  4  shows  the  refractory  telescope. 

Wisconsin  Avenue  passes  just  west  of  the  Observatory  grounds,  run¬ 
ning  approximately  north-south;  see  again  Figure  1.  It  is  the  develop¬ 
ment  along  this  street  that  is  the  concern  of  the  Observatory,  in  par¬ 
ticular,  the  portion  of  Wisconsin  Avenue  that  runs  along  the  northwest 
perimeter  of  the  Observatory.  Development  to  the  southwest  is  already 
proceedi ng. 

There  are  no  buildings  on  the  western  side  of  the  Observatory 
grounds.  Wisconsin  Avenue  and  the  telescopes  are  at  about  the  same  al¬ 
titude.  In  between,  the  terrain  drops  down  some  15  meters  (40  or  50 
feet)  forming  a  little  valley.  The  buildings  housing  the  telescopes  sit 
on  a  small  hill.  The  western  slope  of  this  hill  is  thickly  covered  with 
trees  (Figure  5).  Foliage  on  the  rest  of  the  grounds  is  moderate  (see 
Figures  6  and  7). 
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Figure  2.  Inside  of  T ransit  House  at  the  Naval  Observatory. 


Figure  3.  Outside  the  Transit  House. 
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Figure  4.  The  26-inch  refracting  telescope  used  at  the 
Naval  Observatory. 


Figure  5.  Looking  back  towards  the  woods  covering  the  western  slope 
of  hill  where  the  telescopes  are  located. 


Figure  7.  Looking  from  Wisconsin  Avenue  toward  the  T ransit  House. 


CLIMATE  OF  THE  SITE 


Ambient  air  temperatures  during  a  Washington,  DC  day  in  January  and 
July  were  acquired  from  References  1  and  2.  They  are  plotted  along  with 
the  wall  temperatures  on  many  of  the  figures  included  with  this  report. 

In  January,  the  mean  resultant  surface  wind  is  from  the  northwest 
at  7  mph;  in  July  it  is  from  the  southwest  at  4  mph.  The  wind  direction 
and  velocity  will  be  used  as  parameters  in  this  study.  Construction  of 
concern  to  the  Naval  Observatory  is  to  their  west.  These  analyses  con¬ 
centrate  on  winds  approaching  from  that  general  direction.  Two  wind 
velocities  were  examined:  2  m/sec  (4.5  mph)  and  5  m/sec  (11  mph). 


BUILDING  HEAT  REJECTION 


A  building  can  heat  the  air  in  two  ways,  discharge  from  the  HVAC 
and  convection  heat  transfer  off  the  walls.  HVAC  discharge,  air  condi¬ 
tioning  or  heating,  depends  upon  the  structure  of  the  building,  its 
function,  people  inside,  and  the  weather.  Calculating  HVAC  loads  is 
straightforward,  albeit  tedious.  It  involves  bookkeeping  all  the  above 
variables  and  conducting  mass  and  energy  balances  into  and  throughout 
the  building.  Several  personal  computer  programs  are  available  to  do 
this  bookkeeping.  NCEL  used  the  Elite  Software  Commercial  HVAC  energy 
estimation  and  load  calculation  program  (Ref  3). 

The  role  that  convection  plays  in  heating  the  air  is  best  explained 
by  example.  Figure  8  is  a  schematic  showing  a  cross  section  of  a  build¬ 
ing  wall.  The  arrows  represent  the  different  forms  of  energy  arriving 
at  and  leaving  the  wall:  solar  radiation,  radiation  from  the  wall  to 
the  surroundings,  conduction  through  the  wall,  and  convection  from  the 
wall  to  the  air.  (Conduction  and  convection  can  be  in  either  direc¬ 
tion.)  Applying  conservation  of  energy  to  the  outer  section  of  the  wall 
(brick) , 

Solar  radiation  +  Conduction  from  inside  the  building  -  Radiation 
to  the  surroundings  -  Convection  to  the  air  =  Increase  in  the 
internal  thermal  energy  of  the  wall 

Note  that  this  relationship  is  transient.  Solar  radiation  is  changing 
with  time.  The  temperature  of  the  air  is  also  changing  with  time; 
therefore,  convection  and  radiation  off  the  wall  are  continually  chang¬ 
ing.  Note  also  that  it  is  one  dimensional.  Heat  conduction  along  the 
walls  is  not  included. 

Conservation  of  energy  of  the  wall  can  bp  expressed  in  mathematical 
form.  For  the  brick  (see  Figure  8): 
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e  = 
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of  the  brick,  and 

t  =  time 

x  =  length  variable 
h  =  convective  film  coefficient 


^solar 

a  = 
T  = 

oo 

T.  = 
ins 


incident  solar  radiation 

Stefan-Bol tzmann  constant 
air  temperature 

temperature  of  wall  insulation 


with  equivalent  relationships  for  the  insulation  and  the  inside  of  the 
wall,  resulting  in  four  coupled,  non-linear  partial  c.'i fferential  equa¬ 
tions.  Equations  (1)  were  solved  numerically,  using  a  finite  difference 
program  developed  at  the  NCEL  and  described  in  Appendix  B.  The  program 
includes  a  geometric  analysis  that  takes  the  angle  of  the  sun's  rays  and 
the  shading  imposed  by  the  angle  of  the  rays  into  consideration.  Finite 
subdivisions  of  the  vail,  used  for  these  analyses,  are  shown  on  Figure 
8.  The  temperature  of  the  outer  surface  of  the  brick,  designated  T  , 
is  the  wall  temperature  referred  to  in  this  report. 

Figure  9  was  developed  using  this  program.  It  shows  the  hourly 
variation  of  temperature  of  the  wall  illustrated  on  Figure  8,  assuming  a 
January  day  in  Washington,  DC.  Early  in  the  morning,  air  and  wall  tem¬ 
peratures  are  approximately  equal.  As  the  sun  rises,  the  building  begins 
to  receive  solar  radiation,  and  the  outer  surface  of  the  walls  start  to 
heat  up.  As  the  inside  of  the  building  is  heated,  energy  is  conducted 
through  the  wall  insulation  toward  the  outer  surface,  also  tending  to 
increase  temperature. 

The  wall  temperature  is  now  slightly  higher  than  the  temperature  of 
its  surroundings.  There  is  some  convection  from  the  walls  to  the  air, 
but  of  much  smaller  magnitude  than  the  incoming  solar  radiation,  and  the 
wall  temperature  continues  to  increase.  Heat  is  still  being  conducted 
through  the  insulation  toward  the  outer  part  of  the  walls.  Radiation 
from  the  walls  to  its  surroundings  is  negligible*. 


*Building  walls  never  get  hot  enough  to  re-radiate  a  significant  propor¬ 
tion  of  the  absorbed  solar  radiation. 
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rigure  9.  Building  wall  temperatures  for  a  January  day  in  Washington,  DC. 


Sometime  in  the  middle  of  the  afternoon,  the  wall  temperatures 
reach  a  maximum  and  begin  to  decrease.  Convection  off  the  walls  is  now 
exceeding  the  absorbed  solar  radiation  plus  whatever  heat  is  being  con¬ 
ducted  out  of  the  building.  The  air  temperature  starts  to  decrease  at 
about  this  same  time  accelerating  convection  heat  transfer  off  the 
walls. 

As  night  approaches,  solar  radiation  ceases,  but  heat  convection 
continues.  The  wall  temperature  continues  to  fall,  retarded  slightly  by 
the  small  amount  of  heat  conducted  out  of  the  building.  Finally,  the 
building's  heat  is  shut  off,  only  convection  remains,  and  the  wall  tem¬ 
perature  decreases  on  down  toward  ambient.  The  next  day,  the  cycle  is 
repeated. 

The  situation  is  analogous  in  the  summer  with  one  difference.  The 
inside  of  the  building  would  normally  be  cooler  than  the  outside  part  of 
the  wall.  Heat  would  thus  be  conducted  into  the  building. 

It  is  the  energy  convected  from  the  walls  to  the  air  that  is  of 
interest  here.  Convection  is  represented  mathematically  in  Equations 
(1)  by  the  second  term  in  the  boundary  condition  at  x  =  0, 


Convection  heat  transfer  =  h  (T  .  .  -  T  ) 

s  Ki  n  °° 


(2) 


Note  that  convection  is  proportional  to  the  difference  between  the  wall 
temperature  and  the  temperature  of  the  air.  Consequently,  the  differ¬ 
ence  between  the  wall  temperature  curves  and  the  ambient  temperature 
curve  on  Figure  9  provides  an  indication  of  the  rate  of  convection  heat 
transfer  to  the  air.  Even  as  late  as  midnight,  the  walls  are  heating 
the  air.  The  other  factor  is  the  film  coefficient,  h.  This  factor  var¬ 
ies  with  the  properties  of  the  air  and  its  velocity.  Since  air  proper¬ 
ties  remain  nearly  constant, 


h  ~  Vn 


For  most  wind  velocities  encountered,  n  =s  0.8.  Experimentally  fitting 
the  constant  of  proportionality  for  air  (Ref  2), 


h  =  0.5 


v0.78 


(3) 


Variables  Affecting  Convection  Off  Building  Walls 

A  study  of  the  heat  convected  off  the  walls  of  buildings,  there¬ 
fore,  is  a  study  of  the  variables  affecting  wall  temperature  and  the 
variables  affecting  the  film  coefficient.  In  order  to  make  such  a 
study,  a  baseline  configuration  was  established.  This  configuration  is 
defined  by  Figures  8  and  10  and  Tables  1  and  2  for  the  walls  and  roof. 
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Table  1.  Properties  of  a  Wall 


Material 

Thickness 

(cm) 

Density 
( kg/m3) 

r~  ■  - 

Specific 

Heat 

(J/kg-°C) 

Thermal 
Conductivity 
(W/M  -  °C) 

Brick 

10.0 

1922 

840 

0.7 

Insulation 

10.0 

91 

840 

— 

Bead 

Board 

0.75 

800 

840 

— 

Gypsum 

Board 

0.5 

800 

1090 

0.4 

Table  2.  Properties  of  a  Roof 


Material 

Thickness 

(cm) 

Densi ty 
(kg/m3) 

Specific 

Heat 

(J/kg-°C) 

Thermal 
Conductivity 
(W/M  -  °C) 

Roof i ng 

0.9 

1300 

— 

Concrete 

30.5 

2080 

880 

0.9 

Gypsum 

Board 

1.6 

800 

1090 

0.4 

There  are  two  general  variables:  the  construction  of  the  wall  and 
the  climate/weather.  Three  outer  wall  construction  configurations  were 
examined.  They  are  classified  as  heavy,  medium,  and  light  according  to 
their  heat  capacity  (density  x  specific  heat  x  thickness).  The  baseline 
wall  (brick)  is  classified  as  heavy.  Table  3  gives  examples  of  the  three 
classifications. 

Figure  11  shows  the  west  wall  temperature  variation  during  a  July 
day  with  outer  wall  construction  as  a  parameter.  As  would  be  expected, 
the  light  wall  heats  up  and  cools  off  rapidly.  Convection  off  this  type 
of  wall  would  be  very  high  in  the  afternoon,  low  In  the  evening  and  neg¬ 
ligible  at  night.  Observe  the  change  in  the  difference  between  the  wall 
temperature  and  air  (ambient)  temperature.  The  high  heat  capacity  brick 
wall  heats  up  slowly,  but  there  is  still  convection  off  this  wall  far 
into  the  night. 
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Table  3.  Wall  Classification 


Class 

Example 

Density 
( kg/fti  ) 

Thickness 

(cm) 

Speci f ic 
Heat 

(J/kg-°C) 

Heat 

Capacity*  - 
( J/m  -°Cxl0  q) 

Heavy 

10  cm  brick 

1922 

10 

840 

16.4 

Medium 

Lightweight 

concrete 

961 

10 

840 

8. 1 

Light 

Metal  siding 

7690 

0.6 

420 

2.0 

*Density  x  Specific  Heat  x  Thickness 


In  Figure  11  it  is  assumed  that  all  three  walls  absorb  the  same 
amount  of  solar  radiation.  This  is  not  normally  true.  A  brick  wall, 
for  example,  might  absorb  twice  as  much  of  the  incident  solar  radiation 
as  a  wall  with  metal  siding.  Absorptivity  is  a  measure  of  this;  it  is 
defined  as  the  fraction  of  radiation  striking  a  surface  that  is  absorbed 
by  the  material.  Typical  values  are  shown  on  Table  4.  It  follows  that 
a  wall  with  a  high  absorptivity  will  heat  up  faster.  Figure  12  illu¬ 
strates  the  effect  of  absorptivity  on  building  wall  temperature. 


Table  4.  Values  of  Absorptivity  at 
Solar  Wavelengths  (Ref  10) 


Material 

Absorptivity 

Asphalt 

0.93 

Roofing  tile,  brown 

0.91 

Slate  composition 

0.90 

Red  brick,  red  paint 

0.70 

Concrete 

0.60 

White  enameled  steel 

0.50 

White  marble 

0.45 

Yel low  paint 

0.35 

Aluminum,  polished  plate 

0.29 

Plaster 

0.28 
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Figure  1 1 .  Effect  of  the  construction  on  the  west  wall  temperature  on  a  Washington,  DC 
building  in  July. 


When  the  wall  is  hotter  than  the  inside  of  the  building,  heat  is 
conducted  toward  the  inside  of  the  wall.  The  outer  surface  of  the  wall, 
therefore,  cools  off.  Note  the  location  of  the  insulation  on  Figure  8. 
Changing  the  resistance  of  this  insulation  should  change  the  outer  wall 
temperature.  It  does,  as  shown  on  Figure  13,  but  very  little. 

If  density  perturbations  in  the  air  are  used  as  a  measure,  a  "good" 
wall  is  a  cool  wall,  particularly  in  the  evening  and  at  night.  Heat 
convected  to  the  air  would  then  be  minimum.  A  re-examination  of  Figures 
11,  12  and  13  suggests  that  a  good  wall  should  be  of  light  construction 
and  low  absorptivity.  Wall  resistance  seems  to  have  little  effect. 
Figure  14  compares  the  two  extremes,  combining  all  the  "good"  wall  pro¬ 
perties  and  then  all  the  "bad"  properties,  using  the  west  wall  as  a 
vehicle.  Temperatures  of  the  various  sides  and  roof  of  the  buildings 
represented  by  these  two  extremes  are  plotted  separately  on  Figures  15 
and  16. 

Climate  and  weather  are  not  controllable  variables,  but  they, 
nevertheless,  affect  the  rate  that  heat  is  convected  from  the  walls  of  a 
building  to  the  surrounding  air.  The  velocity  of  the  wind  is  used  as  a 
parameter  on  Figure  17.  Recall  from  the  discussions  of  Equations  (2) 
and  (3)  that  the  rate  of  convection  off  a  flat  surface  is  proportional 
to  both  the  wind  velocity  and  the  difference  between  the  wall  and  the 
air  temperature.  Plots  such  as  Figure  9,  therefore,  only  tell  half  the 
story.  Given  the  same  temperatures,  the  energy  convected  to  the  air  on 
a  windy  day  might  be  several  times  the  energy  convected  on  a  calm  day. 

Consider  Figure  17.  The  incident  solar  radiation  and  the  absorp¬ 
tivity  are  kept  constant.  Consequently,  the  energy  absorbed  during  the 
day  is  the  same  for  all  conditions  examined,  i.e.,  regardless  of  the 
wind  velocity.  It  follows  that  the  total  energy  convected  to  the  air 
over  24  hours  is  also  independent  of  the  wind  velocity. 

The  windier  the  morning,  the  slower  the  walls  heat  up.  Much  of  the 
solar  energy  absorbed  by  the  walls  is  being  immediately  convected  off  to 
the  wind.  Less  energy  is  being  stored  in  the  walls,  however,  and  the 
heat  convected  to  air  during  the  rest  of  the  day  and  night  is  small. 
Conversely,  if  the  air  is  calm,  little  energy  is  convected  to  the  wind. 
The  incident  solar  radiation  is  stored  in  the  walls  as  internal  thermal 
energy,  and  the  walls  heat  up.  It  is  eventually  convected  out  to  the 
wind,  but  it  takes  a  long  time,  into  the  evening  and  possibly  far  into 
the  night. 

If  the  evening  and  night  are  the  time  of  greatest  concern,  the 
worst  condition  occurs  when  a  calm  day  is  followed  by  a  windy  evening 
and  night.  Wall  temperatures  build  up  during  the  day,  and  the  stored 
energy  is  convected  off  at  high  rates  during  the  night.  The  curve 
marked  0/2  on  Figure  17  illustrates  this  possibility.  The  day  is  calm 
until  7  pm;  at  this  time,  a  2  m/sec  (4.5  mph)  wind  starts  up  and  contin¬ 
ues  to  blow  throughout  the  night.  The  film  coefficient  is  the  same  as 
for  the  "2  m/sec"  curve  beneath  it,  but  the  temperature  gradients  are 
about  twice  as  great. 
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Temperature 


re  1 3.  Effect  of  the  overall  wall  resistance  (m2  -  °C/kJ)  on  the  west  wall  temperature 
of  a  Washington,  DC  building  in  July. 
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Temperature  (°F) 


gure  15.  Wall  and  roof  temperatures  assuming  the  best  construction  type. 


Temperature 


Figure  16.  Wall  and  roof  temperatures  assuming  worst  construction  type. 
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re  1 7.  Effect  of  wind  velocity  on  the  west  wall  temperature  of  a  Washington,  DC 
building  in  July. 


Figures  18*  through  21  show  the  variation  in  the  temperatures  of 
all  four  walls  over  a  12-hour  period.  The  time  of  the  year  is  used  as  a 
parameter,  winter  and  then  summer.  A  Washington,  DC  environment  is  as¬ 
sumed.  As  would  be  expected,  the  energy  convected  to  the  air  during  a 
summer  day  is  somewhat  greater.  (Observe  that  the  temperature  differ¬ 
ences  are  greater.)  The  solar  radiation  incident  on  the  walls  is  grea¬ 
ter.  What  would  possibly  not  be  expected  is  the  significant  convection 
to  the  air  that  still  exists  during  the  winter  months.  Notice  also 
that,  for  this  environment,  the  west  wall  is  the  hot  wall  during  the 
summer;  the  south  wall  is  the  hot  wall  during  the  winter.  (The  effect 
of  clouds  and  the  shading  provided  by  trees  and  nearby  buildings  would 
be  similar. ) 

Calculated  Values  of  Wall  Heat  Convection 

Equations  (1)  were  solved  for  a  single  building  and  also  for  the 
entire  portion  of  Wisconsin  Avenue  adjacent  to  the  Naval  Observatory. 
The  dimensions  of  the  buildings  and  their  locations  are  tabularized  in 
Appendix  C.  In  order  to  prevent  the  analyses  from  becoming  overwhelm¬ 
ing,  Wisconsin  Avenue  was  divided  into  eight  regions.  The  HVAC  loads 
and  the  wall  convection  off  buildings  located  in  each  region  were  added 
together,  and  spread  evenly  along  that  portion  of  Wisconsin  Avenue  fall¬ 
ing  within  each  region.  The  concept  is  shown  schematically  on  Figure 
22.  Observe  that  most  of  the  current  development  is  to  the  southwest  of 
the  Observatory. 

Results  for  a  late  afternoon  in  July  are  summarized  on  Table  5. 
Q(cool)  is  the  HVAC  load.  Q(roof)  and  Q(  sides)  are  the  energy  being 
convected  off  the  roof  and  four  sides  of  the  building(s).  Locations 
correspond  to  the  locations  defined  on  Figure  22. 

Discussion  of  Building  Heat  Transfer 

Two  further  observations  need  to  be  emphasized.  First,  convection 
appears  to  be  the  major  source  of  heat  rejected  off  buildings  during 
most  of  the  day  and,  at  least,  into  the  early  evening.  Table  5  compares 
HVAC  and  convection  during  the  late  afternoon.  For  the  single  building, 
convection  is  about  seven  times  greater  than  the  energy  discharged  with 
the  air  conditioning.  For  the  block  sections  of  Wisconsin  Avenue,  the 
difference  is  less,  but  many  of  these  buildings  are  small  one-  and  two- 
story  structures. 

The  second  observation  is  that  convection  (and  HVAC)  off  building 
walls  occurs  into  the  night.  For  buildings  with  walls  of  heavy  construc¬ 
tion,  convection  may  still  be  significant  at  midnight  or  later. 


*Figure  18  is  analogous  to  Figure  9  with  the  absorptivity  increased 
from  0.7  to  0.9. 
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Temperature  (°C) 


NOTES:  (1)  Wind  velocity  =  2  m/sec  (4.5  mph) 

(2)  Absorptivity  =  0.9 

(3)  Overall  wall  resistance  =  0.6  m2-  °C/kJ 


ambient 


Figure  19.  Building  wall  temperatures  in  Washington,  DC  during  a  calm  day  in  July. 
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NOTES:  (1)  Wind  velocity  =  5  m/sec 

(2)  Absorptivity  =  0.9 

(3)  Overall  wall  resistance  =  0.6  ms-  °C/kJ 
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Figure  20.  Building  wall  temperatures  in  Washington,  DC  during  a  windy  day 
in  January. 
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NOTES:  (1)  Wind  velocity  =  5  m/sec 

(2)  Absorptivity  =  0.9 

(3)  Overall  wall  resistance  =  0.6  m2-  °C/kJ 
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Figure  21 .  Building  wall  temperatures  in  Washington,  DC  during  a  windy  day 
in  July. 
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Temperature  (°C) 


Figure  22.  Schematic  of  the  Naval  Observatory  showing  the  division  of  Wisconsin 
Avenue  for  the  purpose  of  lumping  building  loads. 


Table  5.  Summary  of  Calculated  Heat  Rejection  Off  the  Buildings 
Along  Wisconsin  Avenue  in  the  Vicinity  of  the  Naval 
Observatory  During  a  Late  Afternoon  in  July 


Location 

Q  (roof) 
(Btu/hr) 

Q  (sides) 
(Btu/hr) 

Q  (cool) 
(Btu/hr) 

Q  (total) 
(Btu/hr) 

Q  (total) 
(watts) 

1 

2.65E7 

9.95E6 

1.33E7 

4.98E7 

1.46E7 

2 

1.83E7 

1.65E7 

1.47E7 

4.95E7 

1.45E7 

3 

9.73E6 

3.25E7 

4 

1.13E7 

1.22E7 

4.95E7 

1.45E7 

5 

7.16E6 

9.88E6 

3.08E7 

9.02E6 

6 

5.90E6 

2.19E7 

6.42E6 

7 

wm 

4.93E6 

7.48E6 

2.74E7 

8.03E6 

8 

1.19E7 

1.34E7 

7.44E6 

3.27E7 

9.58E6 

Typical 

Single 

Building 

6.54E6 

7.50E6 

2.10E6 

1.61E7 

4.73E6 

Accuracy  of  the  Building  Load  Calculations 

Several  assumptions  were  made  in  order  to  solve  Equations  (1): 

1.  No  shading  from  adjacent  buildings  or  from  trees  was  consider¬ 
ed.  Shading  provided  by  the  building  itself  was  considered. 

2.  Wind  approaches  the  buildings  undisturbed,  with  constant 

velocity,  constant  direction,  and  constant  temperature. 

3.  The  film  coefficient  is  a  function  only  of  the  undisturbed 

wind  velocity.  It  is  the  same  for  all  sides  and  the  roof. 

4.  The  shape  of  the  buildings  was  not  considered.  For  purposes 
of  wall  surface  area  calculations,  all  buildings  were  assumed 
to  be  square. 

5.  Heat  conduction  along  the  walls  is  negligible. 

The  error  introduced  with  the  assumptions  was  not  numerically  as¬ 
sessed.  For  the  Wisconsin  Avenue  configuration,  the  assumption  of  no 
shading  is  good.  Most  of  the  taller  buildings  are  not  located  next  to 
each  other.  Furthermore,  during  the  summer,  the  west  wall  Is  the  hot¬ 
test  wall,  and  It  is  located  along  the  street  where  there  is  no  shading. 
The  roofs  are  not  shaded  regardless. 
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As  discussed  previously,  the  prevailing  winds  in  Washington,  DC  are 
from  the  west.  Buildings  to  the  west  of  Wisconsin  Avenue  are  primarily 
single  family  homes  and  small  shops.  It  is  reasonable  to  assume  that 
these  structures  would  provide  little  resistance  to  the  approaching 
wind.  The  assumption  of  an  undisturbed  wind  is,  therefore,  probably 
good. 

The  next  two  assumptions  are  conservative.  The  upwind  wall,  a 
stagnation  wall,  would  have  a  higher  film  coefficient  than  calculated 
(based  on  the  wind  blowing  over  a  vertical  wall).  The  downwind  wall 
would  have  a  lower  film  coefficient.  During  the  summer,  the  upwind  wall 
i s  the  hot  west  wal 1 . 

A  squire  building  would  have  a  smaller  total  wall  surface  area  than 
a  building  with  the  same  floor  area  but  with  any  other  shape.  The  cal¬ 
culated  total  convection  heat  transfer  to  the  atmosphere  will  be 
smaller. 

Finally,  if  conduction  through  the  walls  induced  by  significant 
temperature  gradients  is  small,  as  shown  on  Figure  13,  conduction  along 
the  walls,  with  very  small  temperature  differences,  is  negligible. 


AIR  DENSITIES 

Wind  blowing  over  the  buildings  is  heated  by  the  HVAC  discharge  and 
by  the  energy  convected  off  the  warm  walls  of  the  buildings.  This  caus¬ 
es  a  local  change  in  the  properties  of  the  air.  These  "pockets"  of  per¬ 
turbed  air  are  carried  along  with  the  wind,  rising  due  to  the  buoyancy 
of  the  now  warmer  air.  The  buoyancy  plus  the  turbulence  of  the  air 
causes  mixing,  which  further  perturbs  the  downwind  structure  of  the  at¬ 
mosphere. 

The  problem  here  is  to  determine  that  structure.  The  behavior  of 
the  atmosphere  can  be  simulated  mathematically  with  a  continuum  repre¬ 
sentation  of  conservation  of  mass  (continuity),  momentum,  and  energy, 
along  with  some  appropriate  relationship  of  state. 

Equation  of  Continuity 

*P  +  3(pu)  +  ilpv).  +  3(p>0  =  o  Mai 

3t  3x  3y  3z  1  J 

Conservation  of  Momentum  (Navier-Stokes  Equations) 


iy  +  u  iy  +  v  ly  +  3U 
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with  similar  expressions  for  the  y  and  z  components. 
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Conservation  of  Energy 
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where:  u,v,w 
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velocity  components  in  x,  y,  z  directions,  respectively, 
den  si ty 

acceleration  of  gravity 

pressure 

vi sco si ty 

body  forces 

thermal  conductivity 

specific  heat 

time 

temperature 


and 


=  velocity  vector  =  ut  +  vj+w£ 
§  =  dissipation  function 


five  equations  and  the  six  unknowns:  the  air  velocity  components  u,  v,  w 
and  the  thermodynamic  properties  of  the  atmosphere,  P,  p,  T.  The  final 
relationship  Is  an  equation  of  state, 

f  (P.P.T)  =  0  (5) 

Turbulence  further  complicates  Equations  (4).  The  transport  prop¬ 
erties  are  now  dependent  upon  the  flow  field.  This  dependency  is  ex¬ 
tremely  complex  and  has,  to  date,  defied  a  theoretical  description.  The 
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usual  approach  is  semi-empirical,  assuming  some  form  for  the  relation¬ 
ships  between  fluid  properties  and  flow  field  and  experimentally  fitting 
constants  in  these  relationships  (Ref  4  and  5). 

The  forces  acting  on  the  flow  field  due  to  the  buoyancy  of  the  warm 
pockets  of  air  must  also  be  included  in  the  mathematical  representation 
of  the  problem. 

Boundary  conditions  (BC)  are  relatively  straightforward.  The  up¬ 
wind  BCs  are  variables  to  be  examined,  wind  velocity,  temperature,  and 
wind  direction.  The  mass  and  energy  input  from  the  buildings  have  been 
discussed  in  detail.  The  geometry  is  set,  i.e.,  the  Naval  Observatory 
and  surroundings.  The  only  BC  decisions  to  be  made  prior  to  attacking 
Equations  (4)  and  (5)  are  to  define  the  type  and  extent  of  the  bounda¬ 
ries.  For  example,  how  much  of  the  atmosphere  must  be  considered?  How 
detailed  must  the  buildings  and  the  ground  be?  These  boundary  condi¬ 
tions  were  established  by  trial  and  error.  They  will  be  discussed 
further. 

Equations  (4)  and  (5)  can  be  simplified  somewhat  by  making  several 
good  additional  assumptions.  The  air  is  considered  incompressible. 
With  wind  velocities  of  a  few  meters  per  second,  density  of  the  air  will 
vary  only  with  temperature.  The  air  is  assumed  to  act  as  an  ideal  gas. 
Finally,  the  analyses  are  considered  independent  of  time.  For  the  minute 
or  two  that  it  takes  the  air  to  pass  between  buildings,  the  temperature 
of  the  walls  will  remain  nearly  constant. 


PHOENICS  PROGRAM 

Even  with  these  assumptions,  Equations  (4)  are  too  complex  to  solve 
in  closed  form.  They  were  solved  numerically,  employing  the  PHOENICS 
code.  PHOENICS  is  a  general  purpose  heat  transfer/f 1 uid  flow  computer 
program  developed  by  CHAM,  Ltd.  of  London  (Ref  6).  It  is  an  iterative 
code,  solving  finite  volume  approximations  of  Equations  (4). 

Turbulence  is  simulated  in  PHOENICS  by  using  an  effective  exchange 
coefficient  consisting  of  both  a  molecular  and  turbulent  component.  For 
example,  the  viscosity  would  be  approximated  as: 

y  -  peff  =  y  +  yturb 

where:  y  =  molecular  (laminar)  viscosity 

vturb  =  v1scos1ty  account  for  turbulence 

Various  empirical  relationships  have  been  used  to  calculate  the  effec¬ 
tive  exchange  coefficients.  PHOENICS  employs  the  k-e  model, 

r2 

^turb*  p  e 

The  parameters  k  (turbulence  kinetic  energy)  and  e  (dissipation  rate) 
then  become  additional  unknowns  to  be  solved  for  (Ref  6): 
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t 


where:  u',  v',  w1  =  velocity  fluctuations 

=  appropriate  mixing  length 

Buoyancy  of  the  hot  air  is  included  in  PHOENICS  by  adding  a  buoyant 
force  to  the  equation  for  conservation  of  momentum  in  the  vertical 
direction  (Ref  7) 


X  = 
where:  g 


-  9  P  8  (T.  -  Tref) 

=  acceleration  of  gravity 


3  =  coefficient  of  expansion 


T  =  ambient  (air)  temperature 
T^  *  =  some  reference  temperature 


Once  the  boundary  conditions  are  set,  PHOENICS  requires  that  the 
domain  of  the  problem  be  divided  up  into  a  series  of  "finite  volumes"  as 
illustrated  on  Figure  23.  Properties  such  as  temperature  and  density 
remain  constant  across  each  volume.  Finite  volume  approximations  of 
Equations  (4)  are  then  solved  for  each  volume,  i.e.,  a  series  of  simul¬ 
taneous  equations,  with  boundary  conditions  for  each  volume  set  by  val¬ 
ues  at  the  adjacent  volumes. 

If  these  finite  volumes  are  infinitesimally  small  and  the  domain  of 
the  problem  is  Infinitely  large,  then  the  PHOENICS  solution  to  Equations 
(4)  is  exact.  Such  an  approach  is  both  numerically  and  economically 
impossible.  Compromises  have  to  be  made.  The  best  approach  is  to  start 
with  a  coarse  grid/small  domain  simulation  that  is  simple  and  cheap, 
then  gradually  decrease  the  grid  size  and  increase  the  domain  until  the 
change  in  the  results  is  less  than  the  desired  accuracy.  ^ 

Time  limitations  prevented  this.  The  fi^e  8  x  8  x  8  rn  grid  shown 
in  Figure  23  was  increased  to  25  x  25  x  25  ni  to  assess  the  grid  size. 
There  was  no  noticeable  change  in  the  range  of  the  calculated  properties 
of  the  atmosphere.  Locations  of  the  isopycnics  (lines  of  constant  den¬ 
sity)  shifted  by  perhaps  10  meters  (30  feet)  in  the  region  directly  over 
the  Observatory. 

Based  on  trial  and  error,  the  atmosphere  remained  undisturbed  300 
meters  above  and  200  meters  to  the  sides  of  the  buildings.  Also  based  on 
trial  and  error,  it  was  found  that  the  actual  shape  of  the  buildings  had 
a  negligible  effect  on  the  downwind  structure  once  beyond  100  meters. 
Buildings  were  assumed  to  be  rectangular  in  shape.  Trees  and  the  slope 
of  the  ground  were  found  to  have  little  effect  on  this  problem. 


RESULTS 

The  objective  of  this  work  is  to  determine  the  downwind  perturba¬ 
tions  of  the  atmosphere.  Results  are  presented  in  the  form  of  iso¬ 
pycnics,  lines  of  constant  density,  predicted  in  the  atmosphere  above 
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the  Naval  Observatory.  Figure  24  shows  the  effects  of  the  single  build¬ 
ing  described  in  Appendix  C.  An  east-west  plane  is  shown.  It  is  July 
at  about  5  pm.  A  2  m/sec  wind  is  blowing  directly  from  the  west,  from 
left  to  right  on  this  figure.  The  air  arrives  at  a  temperature  of 
32.8  °C  (which  corresponds  to  a  density  of  1.151  kg/m'3),  is  heated  by 
the  building,  and  proceeds  on  to  the  transit  house.  The  bottom  part  of 
the  figure  is  identical  to  the  top  with  additional  isopycnics. 

For  an  ideal  gas,  the  density  is  inversely  proportional  to  the  tem¬ 
perature.  Pockets  of  low  density  correspond  to  pockets  of  warm  air. 
The  center  of  the  "plume"  on  Figure  24  is  the  warmest.  It  is  rising, 
spreading,  and  cooling  as  it  is  blown  from  the  building  to  the  transit 
house. 

Figure  25  plots  isopycnics  induced  at  10  pm  by  the  same  building  on 
the  same  day  and  in  the  same  wind. 

The  effects  of  the  entire  portion  of  Wisconsin  Avenue  to  the  west 
of  the  Observatory  are  shown  on  Figure  26.  This  is  Wisconsin  Avenue  as 
it  is  today.  The  difference  between  26(a)  and  26(b)  is  the  direction  of 
the  wind.  On  Figure  26(a),  the  wind  is  blowing  through  the  relatively 
heavily  developed  southwest  portion  of  the  avenue.  On  Figure  26(b),  the 
same  wind  is  blowing  through  the  lightly  developed  northwest  section  of 
Wisconsin  Avenue. 

Figure  27  is  identical  to  Figure  26  with  the  "typical  single  build¬ 
ing"  replacing  all  the  smaller  buildings.  See  Table  C-l. 


DISCUSSION 

The  only  difference  between  Figures  24  and  25  is  the  time.  Both 
show  the  structure  of  the  atmosphere  downwind  from  the  same  single 
building.  The  wall  has  cooled  off  somewhat  between  5  pm  and  10  pm.  As 
a  result,  convection  is  less;  heat  transferred  to  the  air  has  decreased 
from  4.73  million  watts  down  to  2.65  million  watts.  Air  over  the  Obser¬ 
vatory  at  5  pm  is  more  perturbed.  Density  gradients  are  larger.  The 
buoyancy  of  the  warm  pockets  of  air  is  greater.  Therefore,  the  plume  is 
rising  a  little  faster.  Similar  characteristics  could  be  shown  for  the 
other  parameters  examined. 

The  difference  between  Figure  26  and  Figure  27  is  the  size  of  the 
buildings.  Replacing  small  buildings  with  large  buildings,  a  condition 
assumed  to  generate  Figure  27,  will  increase  both  the  HVAC  loads  and  the 
convection  surface  area.  A  more  distorted  atmosphere  can  thus  be  ex¬ 
pected.  The  difference  shows  up  better  when  comparing  Figures  26(b)  and 
27(b);  more  small  buildings  are  being  replaced  in  the  path  of  the  wind. 

The  problem  of  heat  Induced  atmospheric  perturbations  is  very  much 
a  three-dimensional  problem.  Even  for  a  single  building,  there  is  mass 
flow  in  all  directions.  On  Figures  24  and  25,  the  warm  air  is  rising  as 
it  is  swept  eastward  by  the  wind.  As  it  rises,  it  is  replaced  with 
cooler  air  flowing  in  from  the  sides,  1 . e . ,  from  the  north  and  south. 

If  there  are  several  buildings;  for  example,  the  entire  section  of 
Wi  sconsin  Avenue  passing  west  of  the  Observatory;  the  three-dimensional 
flows  become  very  complicated.  Different  buildings  are  heating  the  air 
at  different  rates.  Buoyancy  forces  are  different,  and  different  pockets 
of  air  are  rising  at  different  rates.  Instead  of  relatively  symmetric 
plumes  such  as  shown  on  Figures  24  and  25,  isopycnics  become  very  com¬ 
plex,  impossible  to  intuitively  predict,  compare  Figures  24  and  25  with 
Figures  26  and  27. 
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Figure  23a.  Finite  volume  grid  of  atmosphere  above  the  Naval  Observatory  defined 
for  PHOENICS  -  top  view. 
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Figure  23b.  Finite  volume  grid  of  atmosphere  above  the  Naval  Observatory  defined 
for  PHOENICS  -  side  view. 


afternoon  in  July 
ving  from  the  west 
the  building  -  4.73  (10)6  watts 


induced  by  typical  single  building. 


NOTES:  (1)  E-W  plane  for  a  July  evening  at  10:00  pm 

(2)  Wind  a  2  m/sec,  27.2°C  ,  blow  from  the  west 

(3)  Total  heat  flux  out  of  the  building  =  2.65  (1 0)6  watts 


re  25.  Predicted  isopycnics  (kg/m3)  in  the  vicinity  of  the  Naval  Observatory  induced 


NOTES  (1)  July  evening  at  10:00  pm 
(2)  Wind  at  2  m/sec,  27.2°C 

(a)  SW-NE  plane,  wind  blowing  from  the  SW 
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current  Wisconsin  Avenue  buildings. 


i/ind  blowing  from  SW 


Predicted  isopycnics  (kg/m3)  in  the  vicinity  of  the  Naval  Observatory  induced 
by  buildings  along  Wisconsin  Avenue;  includes  future  development  with 
existing  zoning. 


The  three  dimensionality  of  the  problem  is  better  Illustrated  on 
Figure  28.  This  figure  shows  air  flows  in  the  north-south  plane  through 
the  transit  house.  This  plane  is  perpendicular  to  the  undisturbed  dir¬ 
ection  of  the  wind.  Eddies  with  diameters  of  up  to  200  meters  are  being 
formed,  with  wind  velocities  of  1  m/sec  or  greater  (compared  with  an 
undisturbed  wind  velocity  of  2  m/sec).  The  structure  of  the  atmosphere 
in  this  plane  is  shown  on  Figure  29. 


OBSERVATIONAL  ERROR 

Determination  of  the  observational  error  induced  by  building  heat 
generated  density  perturbations  in  the  atmosphere  was  not  a  part  of  this 
work.  It  was  the  catalyst  for  the  effort,  however,  and  a  short  discus¬ 
sion  is  appropriate. 

The  amount  of  the  anomalous  refraction  depends  upon  the  average 
angle  of  tilt  of  the  isopycnics,  the  thickness  of  the  region  within 
which  the  tilts  occur,  and  the  change  of  the  density  within  this  region. 
A  good  approximation  is  given  by  the  simplified  formula: 

Anomalous  Refraction  =  Constant  x  Average  Tilt  Angle  x  Thickness 

For  the  single  building  and  under  the  conditions  assumed  to  produce 
Figure  24,  this  formula  predicts  an  anomalous  refraction  of  about  0.1 
seconds  of  arc  when  looking  straight  up  from  the  Transit  House.  If  all 
of  Wisconsin  Avenue  is  included,  assuming  future  development,  an  anoma¬ 
lous  error  of  0.3  seconds  of  arc  is  predicted. 

This  seems  like  a  very  small  error  until  one  considers  that  the 
formal  error  of  a  star  position  determined  at  the  Naval  Observatory  is 
about  0.04  seconds  of  arc.  Some  defense  systems  currently  in  use  require 
accuracies  in  the  range  of  0.1  to  0.5  seconds  of  arc. 

This  type  of  observational  error  is  discussed  in  more  detail  in 
Appendix  A  and  in  References  8  and  9. 
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Figure  28.  Predicted  wind  vectors  in  the  vicinity  of  the  Naval  Observatory  when  existing 
Wisconsin  Avenue  buildings  are  each  replaced  with  the  typical  building. 


Appendix  A 


Statement  of 
Dr.  Janes  A.  Hughes 
Director,  Astrometry  Department 
U.S.  Naval  Observatory 
on 

Zoning  Case  No.  87-34 
before  the 

District  of  Columbia  Zoning  Commission 
25  February  1988 

The  Navy,  the  Department  of  Defense  In  general  and  certain  space- 
based  activities  require  sophisticated  celestial  reference  systems  for 
navigation  and  positioning.  Such  systems  consist  of  highly  accurate 
star  positions  and  motions  and  provide  coordinate  grids  which  one  uses 
to  specify  locations  and  directions.  Such  a  grid  Is  entirely  similar  to 
lines  of  latitude  and  longitude  on  the  surface  of  the  Earth,  except  In 
this  case  they  are  located  on  the  sky.  The  word,  "navigation,"  Is  used 
here  In  a  general  sense  which  certainly  Includes  classical  navigation, 
but  which  today  Is  primarily  concerned  with  much  more  advanced  applica¬ 
tions  such  as  satellite  navigation,  missile  guidance  systems,  stellar 
tracking  systems,  satellite  station  keeping  and  various  geodetic  appli¬ 
cations  such  as  determining  azimuths  for  orienting  terrestrial  coordi¬ 
nate  systems.  The  accuracies  of  star  positions  required  by  these  vari¬ 
ous  applications  vary  widely.  Classical  sextant  navigation,  for  exam¬ 
ple,  does  not  require  extreme  accuracy,  but  many  modern  applications 
exhaust  the  currently  available  accuracy.  Naturally  we  are  primarily 
concerned  with  the  most  demanding  requirements.  In  addition,  very 
extensive  academic  and  scientific  use  Is  made  of  star  positions  and 
motions.  In  short,  the  same  data  which  permit  the  specification  of  a 
high  precision  azimuth  for  military  or  geodetic  purposes  may  also  be 
utilized  In  a  purely  academic  way  to  determine  the  distance  to  a  galaxy. 
In  both  cases,  the  crucial  Item  Is  the  existence  of  a  fixed,  non-rotat¬ 
ing  celestial  reference  system.  Such  a  system  Is  often  referred  to  as 
an  Inertial  reference  system. 

The  Astrometry  Department,  of  the  LI.S.  Naval  Observatory  Is  charged 
with  determining  and  maintaining  such  celestial  reference  systems.  The 
U.S.  Naval  Observatory  Is  the  sole  Institution  In  the  United  States 
which  does  this  kind  of  work.  This  Is  so  for  both  historical  and  prac¬ 
tical  reasons.  For  example,  the  work  requires  long-term  commitments  of 
major  resources  which  are  often  Inappropriate  for  universities.  A  sin¬ 
gle  observing  program  can  last  a  decade,  with  a  few  additional  years 
required  for  the  complete  mathematical  and  statistical  discussion  of  the 
data.  This  task  forms  a  major  part  of  the  mission  of  the  Observatory, 
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and  Involves  making  observations  with  various  specialized  telescopes 
Including  Transit  Circle  Telescopes.  It  Is  the  threat  to  the  work  of 
the  Washington  Transit  Circle  which  brings  us  here  as  petitioners. 
These  telescopes  are  small,  but  ultra  precise,  state  of  the  art  instru¬ 
ments  used  to  directly  measure  large  angles  on  the  sky  to  exceedingly 
high  accuracy.  The  observational  data  consists  of  angular  zenith  dis¬ 
tances  measured  to  the  North  and  to  the  South,  plus  the  time  at  which  an 
object  crosses  the  local  meridian.  Such  observations  have  been  made 
from  Washington  at  the  U.S.  Naval  Observatory  for  about  one  century  and 
a  data  base,  internationally  acknowledged  to  be  the  world's  best  of  Its 
kind,  has  been  painstakingly  built  up  from  a  series  of  long  term  observ¬ 
ing  programs.  This  observing  effort  continues  to  the  present  time.  The 
current  observing  program  In  Washington  Includes  about  35,000  stars  plus 
the  Sun,  Moon  and  major  and  minor  planets.  High  quality  astronomical 
observations  of  tens  of  thousands  of  stars  continue  to  be  made  in  Wash¬ 
ington.  Contrary  to  usual  astronomical  practice,  observations  are  made 
around  the  clock,  with  daytime  observations  made  of  the  brightest  stars, 
the  planets  Mercury,  Venus  and  Mars,  and  of  course,  the  Sun.  These  ob¬ 
servations  are  necessary  to  derive  an  Inertial  coordinate  system.  At 
the  present  time  however,  large  commercial  development  around  the  Obser¬ 
vatory,  especially  along  Wisconsin  Avenue,  Is  jeopardizing  these  astro¬ 
nomical  data  gathering  efforts. 

Commercial  development,  as  represented  by  large  buildings,  has  two 
deleterious  effects  on  the  observing  environment:  (1)  it  contributes  to 
the  general  ambient  level  of  artificial  lighting,  and  (2)  It  acts  as  a 
source  of  heat  (In  both  the  summer  and  winter),  which  causes  various 
atmospheric  effects  which  adversely  affect  astronomical  observing.  This 
heating  Is  due  to  space  heating  In  the  winter  and  to  air  conditioning  In 
the  summer,  and  Is  also  due  to  the  heat  "soaked  up"  by  buildings  during 
sunny  days  at  any  time  of  the  year.  All  of  these  heat  sources  are  gen¬ 
erally  proportional  to  the  bulk  of  a  building,  although  specific  usage 
might  produce  relatively  minor  variations.  Dr.  Kodres  of  the  Naval  Civ¬ 
il  Engineering  Laboratory  will  discuss  these  effects  as  part  of  his  pre¬ 
sentation  to  you. 

Light  pollution  Is  currently  of  lesser  concern  to  us.  We  have 
dealt  with  that  problem  by  relocating  certain  observing  programs  to  our 
Flagstaff,  Arizona  facility.  That  station  was  specifically  established 
for  programs  using  larger  telescopes  and/or  requiring  a  dark  sky.  The 
programs  now  being  carried  out  in  Washington  are  minimally  affected  by 
the  generally  high  sky  brightness  In  the  area.  The  Naval  Observatory 
has  deployed  Its  telescopes  to  suit  the  needs  of  the  various  observing 
programs. 

Heat  emitted  from  large  buildings,  however,  remains  a  serious  pro¬ 
blem.  Unlike  the  case  with  light  pollution,  we  cannot  simply  relocate 
the  transit  circle  telescope  program  without  gravely  compromising  a 
vital  data  base  compiled  over  many  years.  It  Is  by  comparing  contempo¬ 
rary  observations  with  observations  made  decades  ago,  that  one  deter¬ 
mines  how  the  stars  are  moving.  These  angular  stellar  motions  are 
caused  by  the  stars'  velocities  through  space.  These  velocities  consist 
of  various  components,  for  example,  galactic  rotation  and  star  streaming 
plus  the  stars'  own  peculiar  Individual  velocities.  In  any  event,  a 
knowledge  of  these  motions  Is  essential  If  one  Is  to  correctly  update 
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the  stellar  reference  system,  or  indeed  to  scientifically  study  the  mo¬ 
tions  themselves.  Combining  new  observations  made  from  a  new  site  of 
unknown  characteristics  would  Introduce  systematic  differences  between 
old  and  new  positions  and  result  In  spurious  motions.  It  is  precisely 
the  homogeneity  of  the  Washington  observations  which  makes  them  the  most 
valued  of  their  type  in  the  world.  This  homogeneity  depends  upon  many 
factors.  For  example,  the  seasonal  motions  of  the  massive  piers  support¬ 
ing  the  telescope  and  the  exact  reactions  of  the  so-called  "meridian 
marks"  which  are  used  to  measure  rates  of  change  of  the  azimuth  of  the 
instrument  are  reliably  known  for  the  Washington  installation  at  the 
level  of  0.01  millimeter  (0.0004  Inch).  This  kind  of  prediction  capa¬ 
bility  comes  strictly  from  many  years  of  study  of  the  Washington  site. 

Heat  energy  Injected  Into  the  atmosphere  by  a  building  causes  vari¬ 
ous  effects.  These  Include  turbulence  and  large  scale  changes  in  the 
structure  of  the  atmosphere  as  represented  by  its  temperature  and  its 
density.  Dr.  Kodres  will  discuss  these  effects  in  more  detail.  From 
the  astronomical  point  of  view,  turbulence  at  great  heights  in  the  atmo¬ 
sphere  (~10  km)  can  cause  rapid  variations  in  the  brightness  of  a  star's 
image  in  a  telescope.  This  effect  is  referred  to  as  "scintillation." 
This  kind  of  variation  is  of  minimal  concern  for  the  kind  of  observa¬ 
tions  being  made  at  the  Naval  Observatory.  In  any  event,  it  is  unlikely 
that  the  nearby  buildings  are  solely  responsible  for  the  observed  scin¬ 
tillation.  In  addition,  turbulence  relatively  nearby  a  telescope  (say 
~1  km  or  nearer)  can  cause  an  image  to  rapidly  and  randomly  move  about 
its  average  position  and  to  vary  in  sharpness.  This  effect  is  called, 
by  astronomers,  "astronomical  seeing."  If  we  were  engaged  in  making 
observations  looking  out  over  some  of  these  large,  nearby  buildings, 
then  we  would  no  doubt  be  concerned  with  the  deterioration  in  the  astro¬ 
nomical  seeing  which  they  would  cause.  However,  we  do  not  observe  over 
these  buildings  so  their  effects  upon  seeing  are  also  not  a  concern. 
Recall  that  the  transit  circle  telescope  observes  only  to  the  North  and 
to  the  South.  Allow  me  to  emphasize  therefore,  that  neither  of  these 
turbulence-induced  effects  are  of  primary  concern  in  the  present  matter. 
The  Observatory  Is  very  concerned,  however,  with  another,  perhaps  more 
subtle  effect,  namely,  "anomalous  refraction."  This  effect  is  caused  by 
the  changes  in  the  atmospheric  structure  mentioned  above. 

To  best  understand  anomalous  refraction  it  is  helpful  to  consider 
what  the  situation  would  be  in  a  normal,  undisturbed  atmosphere.  Con¬ 
sider  the  spherical  Earth  surrounded  by  a  spherical  atmosphere.  The 
structure  of  this  atmosphere  can  be  compared  to  the  structure  of  an 
onion.  That  is,  consisting  of  many  individual  layers  wherein  each  layer 
may  be  considered  to  be  homogeneous.  Since  each  layer  supports  the  lay¬ 
er  above  it,  the  pressure  increases  as  one  moves  downward  in  the  atmo¬ 
sphere  (just  as  in  an  ocean),  and  since  the  density  is  proportional  to 
the  pressure,  the  density  of  the  atmosphere  also  increases  as  one  trav¬ 
els  towards  the  surface  of  the  Earth.  This  change  in  density  causes  a 
light  ray  to  change  its  direction  of  travel  as  It  passes  through  the 
Earth's  atmosphere.  A  change  in  the  direction  of  travel  of  a  light  ray 
under  these  circumstances  Is  described  by  saying  that  the  ray  has  been 
"refracted."  This  refraction  occurs  only  when  a  ray  is  passing  through 
matter  which  has  changes  in  its  density.  The  amount  of  bending  depends 
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upon  the  change  In  density  between  layers  and  upon  the  angle  between  the 
ray  and  the  layers.  A  ray  which  enters  a  layer  perpendicularly  is  not 
refracted  even  If  the  density  changes.  In  any  event,  in  the  case  of  a 
normal  atmosphere,  the  layers  are  regular  and  the  path  of  a  ray  may  be 
calculated  with  high  accuracy.  At  the  Observatory  such  calculations  are 
made  for  every  single  observation  and  the  resulting  bending  is  called 
the  "normal  refraction."  The  final  direction  in  which  a  ray  from  a  star 
is  traveling  when  It  gets  to  a  telescope  (or  to  an  eye  for  that  matter), 
is  the  direction  from  which  the  starlight  seems  to  come  and  this  direc¬ 
tion  gives  the  "position"  of  the  star.  Obviously,  this  position  is  not 
the  true  position  and  therefore  one  must  correct  for  the  refraction  ef¬ 
fect.  The  amount  of  the  normal  refraction  varies  with  ground  level  tem¬ 
perature  and  pressure  as  well  as  with  the  position  of  a  celestial  object 
with  respect  to  the  zenith.  In  general,  the  greater  the  distance  from 
the  zenith  the  greater  the  refraction.  Looking  straight  up,  at  the  ze¬ 
nith,  the  normal  refraction  is  zero  since  the  rays  are  then  perpendi¬ 
cular  to  the  layers  of  equal  density.  Incidentally,  a  line  connecting 
points  in  the  atmosphere  which  have  the  same  density  is  called  an 
isopycnic.  This  is  analogous  to  an  isotherm  which  connects  points  hav¬ 
ing  the  same  temperature,  or  to  an  isobar  connecting  points  of  equal 
pressure.  Perhaps  a  more  familiar  analogy  would  be  the  contour  lines 
which  connect  points  of  equal  elevation. 

Now  we  can  consider  what  happens  If  the  atmosphere  is  not  normal. 
To  help  simplify  the  situation,  we  may  assume  that  over  a  small  area  of 
the  Earth's  surface,  the  curved  layers  of  the  normal  atmosphere  may  be 
considered  to  be  plane  and  parallel  to  the  surface.  When  local,  low- 
level  structural  changes  of  the  atmosphere  are  induced  by  heat  sources, 
these  layers  of  equal  density  may  "tilt."  It  Is  precisely  this  tilting 
which  causes  the  anomalous  refraction.  It  is  called  "anomalous"  because 
normally  the  layers  should  be  plane  and  parallel  and  not  tilted.  A  tilt 
represents  an  anomalous  or  abnormal  situation.  In  addition  to  the  tilt 
effect,  heat  sources  may  also  change  the  value  of  the  local  density  of 
the  air.  This  will  also  change  the  situation,  but  for  the  present  this 
possible  additional  effect  will  be  neglected.  In  any  event  we  take  the 
"total  refraction"  to  be  the  sum  of  the  normal  and  anomalous  refraction, 
and  it  Is  by  this  amount  that  an  observation  must  be  corrected.  The 
anomalous  portion  of  the  total  refraction  presents  a  serious  problem  be¬ 
cause  realistically  it  cannot  be  calculated  like  normal  refraction. 
There  are  too  many  variables.  That  is,  the  amount  of  solar  energy  re¬ 
ceived  by  the  buildings  during  the  day,  the  heat  or  air  conditioning 
usage  by  individuals  or  groups,  the  heat  capacity  of  the  various  build¬ 
ings,  the  geometric  configuration  of  the  heat  sources,  and  so  on.  In 
addition,  the  direction  and  the  speed  of  the  wind  are  crucial  parameters 
which,  for  a  given  set  of  heat  sources,  determine  the  amount  and  direc¬ 
tion  of  the  anomalous  refraction.  The  amount  of  the  anomalous  refrac¬ 
tion  depends  upon  the  average  angle  of  tilt  of  the  Isopycnics,  the 
thickness  of  the  region  within  which  the  tilts  occur  and  the  change  of 
the  density  within  this  region.  A  good  approximation  is  given  by  the 
simplified  formula: 

Anomalous  Refraction  =  Constant  x  Average  Tilt  Angle  x  Thickness 
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Compared  to  the  amount  of  normal  refraction,  the  anomalous  refrac¬ 
tion  is  not  large.  However,  compared  to  the  current  observational  accu¬ 
racy  and  the  required  quality  of  the  astronomical  data,  anomalous  re¬ 
fraction  can  easily  be  a  very  significant  contributor  to  the  total 
error. 

At  this  point  it  is  useful  to  consider  the  size  of  the  angular 
quantities  which  are  being  discussed.  In  other  words,  i-r,  the  present 
context,  what  is  "large?"  What  is  the  size  of  the  observational  error? 
What  accuracy  ’ s  required?  We  are  dealing  with  angles  here  and  everyone 
is  familiar  with  the  concept  of  an  angular  degree.  Ninety  degrees  make 
up  what  is  called  a  right  angle.  A  degree  is  subdivided  into  60  minutes 
of  arc,  not  to  be  confused  with  minutes  of  time.  In  turn,  a  minute  of 
arc  is  subdivided  into  60  seconds  of  arc,  so  that  a  degree  contains 
60  x  60  =  3600  seconds  of  arc  and  a  right  angle  contains  324,000  seconds 
of  arc.  These  quantities  are  written  as  e.g.,  ninety  degrees  =  9C°,  35 
minutes  of  arc  =  35',  and  25.13  seconds  of  arc  =  2 5 V 1 3 ,  and  so  on.  To 
visualize  what  a  second  of  arc  is,  consider  a  dime.  In  order  to  have  a 
dime  appear  to  be  one  second  of  arc  in  angular  size,  it  would  have  to  be 
observed  from  1,000  feet  away,  but  edgewise!  That  is  a  second  of  arc. 
The  normal  refraction  at  the  zenith,  as  stated  earlier,  is  zero,  but  45° 
from  the  zenith  it  is  about  58V5.  Near  the  horizon,  the  normal  refrac¬ 
tion  is  about  30'.  Now  the  error  of  a  single  observation  with  the  Tran¬ 
sit  Circle  Telescope  is  about  0 V 2 ,  and  the  dime  must  now  be  placed  one 
mile  away,  still  looking  at  it  edgewise.  By  combining  multiple  obser¬ 
vations  statistically,  the  formal  error  of  a  star  position  determined  in 
this  way  at  the  Naval  Observatory  in  Washington  is  about  0'.'04,  and  the 
dime  is  now  five  miles  away,  still  on  edge. 

Who  cares  about  such  tiny  angles?  Well,  some  defense  systems  cur¬ 
rently  in  use  require  accuracies  in  the  range  of  0V 1  to  0'.'5  for  their 
calibration  accuracy.  The  calibrator,  namely  the  stars,  must  be  signif¬ 
icantly  more  accurate  in  order  to  be  used  as  a  standard.  That  is,  the 
calibrator  must  not  contribute  any  significant  additional  error  to  the 
overall  error  of  the  system.  This  is  why  the  accuracy  now  being  achieved 
is  barely  meeting  requirements.  It  should  be  understood  that  operational 
systems  are  being  discussed,  not  future  schemes  or  drawing  board  extra¬ 
polations. 

From  the  point  of  view  of  astronomical  science,  higher  accuracies 
are  always  needed  to  expand  our  knowledge  of  the  universe.  Even  so,  our 
request  for  protection  has  not  been  generated  to  simply  satisfy  our  in¬ 
tellectual  curiosity  or  even  as  a  response  to  traditional  academic 
desires,  regardless  of  the  very  real  value  of  such  things.  No,  our  need 
for  protection  has  arisen  from  concern  about  very  real  and  practical 
national  requirements  which  can  only  be  satisfied  by  the  Naval  Observa¬ 
tory.  I  remind  you  that  in  the  United  States,  it  is  only  the  Naval 
Observatory  which  makes  fundamental  astrometric  observations.  For  the 
record,  may  I  state  that  the  U.S.  Naval  Observatory  is  proud  of  the  con¬ 
tributions  it  has  made  to  astronomical  science  and  anticipates  with 
pride  the  future  contributions  it  hopes  to  make. 

With  the  foregoing  as  background,  some  typical  effects  of  large 
buildings  will  now  be  exhibited  and  investigated. 
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First,  a  building  is  characterized  by  such  things  as  its  bulk,  the 
area  of  its  exterior  walls  and  roof,  building  materials,  R  values,  HVAC 
energy  usage,  and  so  on.  From  this  data  an  estimate  of  the  wall  and 
roof  temperatures  under  specific  ambient  conditions  is  derived.  These 
ambient  conditions  include  the  outside  and  inside  air  temperatures  and 
the  angle  and  intensity  of  the  incoming  solar  radiation  with  any  shading 
effects.  The  contributions  to  the  total  heat  from  people,  light,  and 
infiltration  are  also  included.  The  computer  program  used  for  the  cal¬ 
culations  is  very  general.  It  calculates  the  amount  of  solar  energy 
incident  upon  each  wall  and  upon  the  roof,  the  amount  of  energy  re¬ 
radiated,  conducted  through  the  walls,  convected  to  the  wind  and  stored 
within  each  wall.  In  this  way  a  realistic,  dynamical  picture  of  the 
building  as  a  heat  source  is  generated. 

Using  the  data  described  above,  the  resulting  air  flow  is  then  cal¬ 
culated  by  a  very  sophisticated  program  which  numerically  solves  the 
applicable  equations  of  hydrodynamics,  the  so-called  Navier-Stokes  equa¬ 
tions.  Dr.  Kodres  will  describe  this  process  and  also  give  additional 
details  regarding  the  characterization  of  a  building. 

Additional  input  data  for  this  program  includes  the  direction  and 
speed  of  the  wind,  the  relative  locations  of  the  buildings  and  the  tele¬ 
scope  and  so  on.  The  output  of  this  program  includes  the  value  of  the 
density  as  it  varies  in  three  dimensions.  From  these  data,  cross  sec¬ 
tions  of  the  isopycnics  are  then  plotted.  It  is  from  these  plots  that 
one  can  calculate  or  estimate  the  amount  of  the  anomalous  refraction  in 
a  particular  direction.  The  entire  process  is  time  consuming  and  com¬ 
plicated.  This  is  because  of  the  scientific  rigor  of  the  physics  used 
and  the  completeness  of  the  mathematical  approach.  The  computations 
were  carried  out  in  part  on  a  CRAY  computer.  This  computer  is  acknow¬ 
ledged  to  be  the  fastest  and  most  sophisticated  machine  available  today. 
In  short,  the  results  which  are  being  presented  to  you  are  products  of 
the  best  technology  available  today.  The  Observatory  recognizes  its 
responsibilities  in  this  serious  matter,  and  it  has  conscientiously  made 
every  effort  to  properly  derive  and  interpret  the  data  presented  to  you 
today.  For  example,  the  use  of  the  full  Navier-Stokes  equations  compli¬ 
cates  the  analysis  greatly,  but  these  are  the  rigorous  equations  which 
describe  the  situation.  Anything  else  is,  by  definition,  an  incomplete 
descri pti on . 

First,  let's  consider  a  single  building.  For  purposes  somewhat 
parallel  to,  but  essentially  independent  of  the  present  hearing,  a  par¬ 
ticular  proposed  Planned  Unit  Development  has  been  investigated  in  some 
detail.  In  this  case  the  plot  depicts  the  isopycnics  as  they  occur  in 
the  East-West  plane.  The  building,  in  this  case,  is  essentially  West  of 
the  Transit  Circle  telescope.  The  wind  is  taken  to  be  from  the  West  at 
2  m/s  (~4  mph).  The  time  is  a  late  afternoon  in  July.  The  tilt  of  the 
isopycnics  is  obvious.  Using  the  formula  given  earlier,  one  calculates 
an  anomalous  refraction  (when  looking  straight  up)  of  about  O'.'l  in  this 
case.  This  is  200%  of  the  standard  error  of  the  Transit  Circle,  0 V 0 4 , 
mentioned  earl ier. 

We  are  not  dealing  with  a  single  building,  however,  and,  therefore, 
we  now  consider  the  effects  of  the  large  buildings  currently  located  on 
Wisconsin  Avenue.  About  32  buildings,  ranging  in  location  from  the  1800 
block  of  Wisconsin  Avenue  up  to  the  intersection  of  Calvert  Street  and 
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Wisconsin  Avenue,  were  included  in  the  calculation.  Here  I  will  show 
the  results  of  the  calculations  made  assuming  a  SW  wind  and  WNW  wind, 
for  the  same  time  as  the  preceding  example.  In  the  case  of  the  SW  wind, 
the  anomalous  refraction  (again  looking  straight  up)  is  estimated,  from 
the  plotted  isopycnics,  to  be  about  0'.'05.  With  the  wind  from  the  WNW 
the  estimated  anomalous  refraction  drops  to  around  0'.'015.  Since,  at  the 
present  time,  the  largest  buildings  are  located  more  to  the  SW  of  the 
telescope,  the  difference  between  these  two  results  is  exactly  what  one 
would  expect,  and  thus  gives  great  confidence  in  the  correctness  of  our 
calculations.  It  is  evident  that  we  already  have  a  problem  when  the 
wind  is  from  the  SW.  Of  course,  if  even  larger  buildings  were  placed  to 
our  SW,  the  problem  would  become  even  worse.  However,  to  illustrate  the 
virtue  of  our  proposed  overlay  district,  consider  the  case  of  a  WNW  wind 
with  additional  development  more  directly  West  and  Northwest  of  the 
telescope.  By  "additional  development"  is  meant  that  which  could  occur 
under  present  zoning,  including  the  PUD  process.  The  anomalous  refrac¬ 
tion  now  rises  to  about  0 V 3 .  On  the  other  hand,  if  one  applies  the 
overlay  concept,  then  the  anomalous  refraction  increases  to  only  about 
O'.1 06 . 

In  order  to  present  our  case  to  you  it  has  been  necessary  to  dis¬ 
cuss  many  concepts  and  effects  which  surely  have  never  been  an  element 
in  a  zoning  petition  before.  Astronomical  observations  are  bad  enough, 
let  alone  something  with  the  unlikely  name  of  "anomalous  refraction," 
not  to  mention  that  verbal  absurdity,  "isopycnic."  We  are  perfectly 
aware  that  such  things  just  don't  form  a  part  of  most  people's  daily 
experience.  Nevertheless,  these  things  are  very  real.  The  rigorous  ap¬ 
proach  which  we  have  used  (in  lieu  of  simple  approximations)  has  given 
results  which  are  very  real.  The  uses  to  which  the  results  of  our  as¬ 
tronomical  observations  are  put  are  very  real.  And  yes,  the  threat  to 
our  continued  ability  to  satisfy  national  requirements  is  also  very 
real . 

Finally,  I  would  like  to  say  that  the  methods  and  principles  used 
by  us  have  been  questioned  by  some.  For  example,  the  use  of  the 
Navier-Stokes  equations  has  been  criticized.  This  criticism  is  untena¬ 
ble  and  frankly  would  be  laughed  out  of  existence  if  submitted  to  any 
fluid  flow  expert  in  the  world.  Apropos  of  this,  we  stand  ready  to  sub¬ 
mit  our  scientific  results  to  any  acknowledged  independent  expert,  and 
incidentally,  Dr.  Kodres  and  I  plan  to  submit  a  paper  on  our  results  to 
an  appropriate  learned  journal  for  judgment  by  our  peers.  It  must  be 
made  known  that  to  date,  none  of  our  disparagers  has  made  a  single  cal¬ 
culation  of  a  single  isopycnic  much  less  the  resulting  refraction.  On 
occasion  some  very  simplified  building  models  have  been  calculated  by 
them  on  a  home,  personal  computer.  Such  calculations  may  suffice  for 
crude,  "back  of  the  envelope"  estimations,  but  we  believe  that  the  Zon¬ 
ing  Commission  and  interested  citizens  deserve  better.  We  have  given 
this  problem  our  best  scientific  "shot,"  and  we  request  your  help  in 
resolving  the  serious  problems  which  our  scrupulous  calculations  have  so 
clearly  demonstrated. 
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Appendix  B 


Computer  Program  to  Simulate  the  Thermal 
Characteristics  of  the  Walls  and  Roof 
of  Buildings 

by 

Cynthia  H.  Ruf 

A  computer  program  was  developed  at  the  Naval  Civil  Engineering 
Laboratory  to  determine  the  temporal  and  spacial  variation  in  building 
wall  and  roof  temperatures.  The  program,  written  in  BASIC,  solves  the 
finite  difference  approximation  of  the  energy  equation  governing  heat 
transfer  to  and  through  the  building  walls. 

The  program  considers  the  amount  of  radiation  incident  upon  the 
building,  taking  into  account  the  angle  of  the  sun's  rays  and  the 
shading  imposed  by  the  angle  of  the  rays.  It  also  considers  the  amount 
of  energy  reradiated,  conducted  through  the  wall,  convected  to  the  wind, 
and  stored  within  the  wall  as  internal  thermal  energy.  The  simulation 
has  the  capability  to  vary  the  thermal  resistance  of  the  wall  in  order 
to  determine  the  sensitivity  of  wall  temperatures  to  wall  construction. 
This  program  also  allows  the  user  to  vary  the  convection  film 
coefficient  in  order  to  examine  the  sensitivity  of  wall  temperatures  to 
wind  speed. 

This  Appendix  is  the  user's  guide  to  the  program.  In  addition,  the 
BASIC  listing  and  typical  output  are  provided. 


USER'S  MANUAL 

A  BASIC  computer  program  was  developed  to  calculate  building  skin 
temperatures  due  to  solar  loads.  This  program  was  specifically 
developed  to  determine  these  temperatures  in  order  to  analyze  the 
thermal  interference  generated  by  a  building.  The  program  was  developed 
for  a  specific  site  and  therefore  is  not  particularly  user  friendly  or 
interactive.  However,  with  the  background  information  presented  here, 
the  user  should  be  able  to  modify  the  necessary  inputs  and  design 
parameters  for  any  new  situation. 

The  program  was  developed  from  segments  of  an  existing  BASIC 
program  which  contained  extensive  solar  calculations.  These 
calculations  were  modified  and  used  in  conjunction  with  additional  heat 
transfer  calculations  to  determine  the  building  skin  temperatures. 

In  general,  the  program  calculates  the  total  solar  load  for  walls 
of  specified  orientation  for  one  representative  day  of  each  month  of  the 
year.  The  program  considers  solar  geometry,  weather,  location,  and 
several  variable  design  parameters  in  its  calculations.  The  variable 
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design  parameters  allow  the  user  to  determine  best/worst  case  scenarios 
for  the  specific  problem  under  study. 

The  program  is  basically  divided  into  three  segments:  input  data, 
solar  insolation  calculations  and  heat  transfer  calculations.  The 
program  uses  equations  from  the  ASHRAE  Applications  Handbook  (Ref  1)  as 
well  as  other  methods  for  its  solar  insolation  calculations.  The 
following  is  a  brief  description  of  the  information  necessary  in  the 
input  data  section. 

The  information  specified  in  lines  100-300  are  the  data  used  to 
calculate  direct  normal  insolation  at  the  earth's  surface  on  a  clear 
day.  'a'  is  the  apparent  solar  irradiation  at  air  mass  zero  for  each 
month;  ' b '  is  the  atmospheric  extinction  coefficient;  and  1  c '  is  the 
ratio  for  the  diffuse  radiation  on  a  horizontal  surface  to  the  direct 
normal  irradiation  (Ref  1). 

Next,  the  julian  date  for  one  representative  day  of  each  month  is 
specified  in  DAY,  line  330. 

PAS,  line  390,  is  the  percent  of  available  sunshine  at  the 
location. 

The  data  for  KTBAR,  line  450,  represents  a  measure  of  cloudiness 
and  other  atmospheric  conditions  which  attenuate  solar  radiation  at  a 
given  location.  Site  specific  values  for  KTBAR  can  be  obtained  from 
Reference  11. 

The  monthly  average  ambient  temperature  is  specified  in  TAVG,  line 
510.  The  percentage  of  the  daily  temperature  range  is  specified  in  line 
570,  PCENT.  PCENT  is  used  along  with  DRANGE,  line  810,  which  is  the 
daily  temperature  range,  to  determine  the  ambient  temperature  throughout 
the  day  for  the  location. 

The  TILT  parameter,  line  630,  must  be  set  to  90  degrees  for  a 
perpendicular  wall  or  to  0  degrees  for  a  horizontal  surface  such  as  a 
roof.  TILT  can  also  be  set  to  any  intermediate  value  if  needed. 

The  latitude  (LAT),  longitude  (LONG),  and  time  zone  (TZN)  for  the 
location  are  specified  in  lines  680-700. 

In  addition,  to  specify  the  tilt  of  the  surfaces,  the  user  must 
also  specify  the  orientation,  i.e.,  N,  S,  SW,  etc.  A  subroutine  which 
calculates  the  solar  angles  for  each  orientation  is  also  included  in  the 
program,  lines  4000-4410. 

Some  other  design  parameters  included  in  the  program  are  PHO,  the 
ground  reflectance;  TIN,  the  indoor  design  temperature;  and  HCONV,  the 
convective  heat  transfer  coefficient. 

The  wall  model  that  is  used  in  the  program  consists  of  layers  of 
materials  with  variable  resistances,  including  resistances  for 
convective  losses  at  both  the  inside  and  outside  wall  interfaces.  The 
design  densities,  thicknesses,  and  specific  heats  of  the  materials 
incorporated  in  the  wall  model  are  included  in  the  design  parameters, 
lines  990-1070.  Similarly,  the  appropriate  resistances  for  the  model 
wall  layers  are  specified  in  lines  1080-1210. 

Once  all  of  the  necessary  information  for  the  area  to  be  studied  is 
entered  in  the  program,  the  program  performs  some  preliminary 
calculation.  The  solar  declination  is  calculated.  Sunrise  and  sunset 
times  are  determined  in  line  1520  for  use  within  the  program.  The  daily 
totals  for  convective,  conductive,  solar,  and  radiative  heat  are 
initialized  for  each  new  day,  lines  1580-1680.  The  ambient  temperature 
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is  calculated  in  line  1700  from  the  average  ambient  temperature,  the 
daily  range,  and  the  wall  outside  surface  temperature  (TSKIN)  to  the 
ambient  temperature. 

The  program  then  performs  calculations,  using  15-minute  increments, 
TSTEP  =  0.25,  over  the  course  of  the  specified  day.  The  program  calcu¬ 
lates  various  solar  angles,  such  as  the  azimuth  and  the  solar  incident 
angle,  which  are  used  to  determine  solar  insolation.  These  angles, 
along  with  the  previously  describes  input  data,  are  then  used  to  calcu¬ 
late  the  insolation  on  the  building  walls.  The  insolation  calculations 
begin  on  line  2040. 

Next,  the  program  performs  the  heat  transfer  calculations  beginning 
on  line  2300.  Conduction,  convection,  and  radiation  are  all  included  in 
the  analysis.  The  conductive  heat  transfer  is  determined  in  subroutine 
5000,  where  the  program  determines  the  energy  stored  in  each  layer  of 
the  wall  model.  The  convective  heat  transfer  is  determined  by  the  skin 
temperature,  the  ambient  temperature,  and  the  design  convective  heat 
transfer  coefficient,  line  2330.  The  radiation  transferred  into  the 
wall,  calculated  above,  is  the  total  solar  insolation.  The  radiation 
transferred  out  from  the  wall  is  assumed  to  radiate  to  surrounding 
buildings  with  temperatures  similar  to  those  of  the  wall  under  study, 
and  is  specified  in  line  2350.  The  wall  skin  temperature  can  then  be 
determined  in  subroutine  6000  where  the  program  determines  the 
temperature  at  each  model  wall  layer  interface. 

The  program  outputs  the  wall  temperature,  in  15-minute  increments 
throughout  the  day.  The  program  will  then  loop  to  the  next  month.  The 
program  can  also  be  modified  to  loop  to  the  next  wall  orientation  to 
complete  the  four  (or  more)  walls  of  the  structure,  as  well  as  the  roof. 

A  sample  output  for  the  following  BASIC  program  is  included. 


PROGRAM  LISTING 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

391 

180 

190 

200 


rem  This  is  a  basic  program  to  calculate  the  skin 
rem  temperature  of  a  building  due  to  solar  loads 
rem 
rem 
rem 

dim  a(12),  b(12),  c(12),  day(12),  drange(12) 

dim  pas(12),  ktb(12),  tavg(12),  pcent(24) 
rem  *************************  +  ***-,********  +  ***+■*** 

Input  Data  Section 

★  ★a***********************'*'-*'*'*-***'**-*'****'*'*'**** 


Dimension  statements 


rem 

rem 

rem 

rem 

rem 

rem 

rem 

rem 


a  =  the  apparent  solar  irradiation 
b  =  the  atmospheric  extinction  coefficient 
c  =  the  ratio  of  the  diffuse  radiation  on  a 
horizontal  surface  to  the  direct  normal 
irradiation  (from  ASMRAE) 


data  390,  385,  376,  360,  350,  345,  344,  351.  365,  378,  387, 


for  q  =  1  to  12 
read  a(q) 
next  q 


210  rem 

220  data  0.142,  0.144,  0.156,  0.180,  0.196,  0.205,  0.207,  0.201, 
0.177,  0.160,  0.149,  0.142 
230  for  g=l  to  12 
240  read  b(g) 

250  next  g 
260  rem 

270  data  0.058,  0.060,  0.071,  0.097,  0.121,  0.134,  0.136,  0.122, 
0.092,  0.073,  0.063,  0.057 
280  for  z  =  1  to  12 
290  read  c(z) 

300  next  z 
310  rem 

320  rem  **********  day  0f  the  year  ********** 

330  data  16,  47,  75,  106,  136,  167,  197,  228,  259,  289,  320, 

350 

340  for  d  =  1  to  12 
350  read  day(d) 

360  next  d 
370  rem 

380  rem  **********  percent  available  sunshine  ********** 

390  data  0.75,  0.73,  0.84,  0.86,  0.89,  0.93,  0.96,  0.95,  0.92, 

0.91,  0.85,  0.80 
400  for  W  =  1  TO  12 
410  read  pas(w) 

420  next  w 
430  rem 

440  rem  **********  Ktbar  ********** 

450  data  0.404,  0.435,  0.450,  0.47,  0.488,  0.514,  0.503,  0.494, 
0.489,  0.475,  0.417,  0.374 
460  for  n  =  1  to  12 
470  read  ktb(n) 

480  next  n 
490  rem 

500  rem  **********  ave  month  ambient  temps  ********** 

510  data  32.1,  33.8,  41.8,  53.1,  62.6,  71.1,  75.3,  73.6,  66.9, 
55.9,  44.7,  34.0 
520  for  x  =  1  to  12 
530  read  tavg(x) 

540  next  x 
550  rem 

560  rem  **********  percent  of  the  dailv  temp  range  ********* 
570  data  0.87,  0.92,  0.96,  0.99,  1.00,  0.98,  0.93,  0.84,  0.71, 
580  data  0.58,  0.68,  0.76,  0.82 
590  for  1  =  1  to  24 
600  read  pcent(l) 

610  next  1 
620  rem 

630  rem  **********  other  data  ********** 

640  rem  the  tilt  must  be  90  for  a  wall 
650  rem  the  tilt  must  be  0  for  a  roof 
660  tilt  =  90 
670  rem 


B-4 


********* 


680  1  a t  =  38 

690  long  =  77 
700  tzn  =  5 
710  rem 

720  rem  *********  direction  of  the  wall 
730  rem  direction  =  0  for  north 

740  rem  direction  =  1  for  east 

750  rem  direction  =  2  for  west 

760  rem  direction  =  3  for  south 

770  rem  direction  =  4  for  roof 

790  rem 

800  rem  **********  daily  temperature  range  ********** 

810  data  24.2,  22.7,  22.7,  21.7,  21.0,  19.8,  18.8,  19.0,  22.6, 

25.8,  28.6,  25.2 
820  for  i  =  1  to  12 
830  read  drange(i) 

840  next  i 

850  rem 

860  rem  **********  ground  reflectance  ********** 

870  rho  =  0 
880  rem 

890  rem  **********  design  parameters  ************ 

900  rem  ***********  design  indoor  temperature  *********** 

910  tin  =  75 
920  rem 

930  rem  ************  design  convective  heat  transfer  coef  ***** 

940  hconv  =2.36 
950  rem 

960  rem  ***********  design  emissivity  of  outer  surface  ************ 
970  eout  =  0.75 

980  rem 

990  rem  ******  design  density  of  model  wall  surfaces  ***** 

1000  dll  =  120 

1001  dl2  =  5.7 

1002  d 1 3  =  50 

1003  dl4  =  50 
1010  rem 

1020  rem  *****  design  thickness  of  model  wall  surfaces  ***** 

1030  111  =  0.83 

1031  112  =  0.83 

1032  113  =  0.0625 

1033  114  =  0.052 
1040  rem 

1050  rem  *****  design  specific  heat  of  model  wall  surface  ***** 

1060  shl 1  =  0.2 

1061  shl2  =  0.2 

1062  shl3  =  0.2 

1063  sh 14  =  0.26 
1070  rem 

1080  rem  *********  the  r  -  values  for  each  wall  layer  **** 

1090  rem  ***  r  -  values  can  be  found  in  ASHRAE  HOF  chapter  26 
1100  rem 

1110  rlayer  =  0.11 
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1120  rlayer  =  4.25 
1130  rl  =  3.9 
1140  r2  =  0.56 
1150  rin  =  0.685 
1160  layerl  =  4 
1170  layer2  =  4 

1200  rem  ************************************************ 

1210  rtot  =  ( layerl*r1ayerl)  +  (Iayer2*rlayer2)  +  ri  +  r2 
1220  rem  ************************************************ 

1230  rem  *********  the  Stefan-Bol tzmann  constant  ****** 

1240  sigma  =  1.74E-09 
1250  rem 

1260  lat  =  lat  *  0.01745 
1270  long  =  long  *  0.01745 
1280  tilt  =  tilt  *  0.01745 
1290  stilt  =  si n( tilt) 

1300  ctilt  =  cos(tilt) 

1310  slat  =  sin(lat) 

1320  clat  =  cos(lat) 

1330  rem 

1340  rem  *********  set  the  time  increment  ********* 

1350  tstep  =  0.25 

1360  rem  ********************************************* 

1365  for  direction  =  0  to  3 

1370  rem  ********************************************* 

1380  rem  **  begin  the  monthly  calculations  ** 

1390  rem 

1400  rem  for  mon  =  1  to  12 

1410  for  mon  =  1  to  8  step  6 

1411  rem 

1420  decl  =  23. 45* si n(360*. 01745 *(284+day( mon)) /365)*. 01745 
1430  sdecl  =  sin(decl) 

1440  cdecl  =  cos(decl) 

1450  rem 

1460  rem  **********  column  headers  ***************** 

1470  print  "month",  "day  o  the  year" 

1480  print  mon,  day(mon) 

1490  if  (direction  =  0)  then  print  "the  wall  faces  north" 

else  if  (direction  =  1)  then  print  "the  wall  faces  east" 

else  if  (direction  =  2)  then  print  "the  wall  faces  west" 

else  if  (direction  =  3)  then  print  "the  wall  faces  south" 

else  if  (direction  =  4)  then  print  "this  is  the  roof" 

1500  print 
1510  print 

1520  rem  ********calculate  sunrise  and  sunset  ********** 

1530  x  =  (-tan(lat)*tan(decl )) 

1540  sunris  =  -atn(x/sqr(-x*x+l))+3. 1416/2 
1550  irise  =  13-( sunris/. 2618) 

1560  i set  =  12+( sunris/. 2618) 

1570  rem 

1580  rem  ********  initialize  the  daily  totals  **************** 
1590  qday  =  0 
1600  qaday  =  0 
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1610  Ibttot  =  0 
1620  idttot  =  0 
1630  itttot  =  0 
1640  dayhrs  =  0 
1650  tqradout  =  0 
1660  tqcond  =  0 
1670  tqconv  =  0 
1680  rem 

1690  rem  *******  calculate  the  ambient  temperature  ******8 

1700  tamb  =  (tavg(mon)+  drange(mon)/2)  -  pcent  ( i ri se)*drange(mon) 

1701  tskin  =  tamb 

1702  tl  =  tamb 

1703  t2  =  tamb 

1704  t3  =  tamb 

1705  t4  =  tamb 

1706  t5  =  tamb 

1707  t6  =  tamb 

1708  t7  =  tamb 

1709  t8  =  tamb 

1710  t9  =  tin 

1711  tlO  =  tin 

1712  tmax  =  tskin 

1720  rem  ********  set  the  time  limits  ***************** 

1730  i  =  int(irise/tstep) 

1740  n  =  int(iset/tstep) 

1750  rem 

1760  rem  *********  output  headers  ***************** 

1770  print  "time",  "tamb",  "tskin" 

1780  rem 

1790  rem  ******************************************* 

1770  print  "time",  "tamb",  "tskin" 

1780  rem 

1790  ************************************************ 

1800  rem  begin  the  daily  calculations 

1810  for  j  =  i  to  n 

1820  k  =  j*tstep 

1830  xi  =  ( 12-k)* . 2618 

1840  rem 

1850  rem  *■■**“•''*"  solar  angles  ■**„*•* ********** 

1860  sbeta  =  (clat*cdecl*cos(xi)+slat*sdecl ) 

1870  if  (sbeta  <  .017)  then  sbeta  =  .017 
1880  cbeta  =  sqr( l-sbeta*sbeta) 

1890  if  (cbeta  =  0)  then  obeta  =  .001 
1900  rem 

1910  rem  ********  solar  azimuth,  psi  *************** 

1920  xpsl  =  (cdecl*sin(x1 )/cbeta) 

1930  psi  =  atn(xpsi/sqr(-xpsi*xpsi+l)) 

1940  cosh  =  cos( . 2618*(k-12+tzn)-long) 

1950  test  =  sdecl*clat/(cdecl*s1at) 

1960  if  ((cosh<=test)  and  ( k< 12 ) )  then  psi  =  3.1416  -  psi 
1970  if  ((cosh<=test)  and  (k>12))  then  psi  =  -3.1416-psi 
1980  gasub  4000 
1990  cazm  =  co5(azm) 
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2000  ctheta  =  cbeta*cazm*sti It  +  sbeta*ctilt 
2010  if  (ctheta  <=  0)  then  ctheta  =  0 

2020  rem  ************************************************** 

2030  rem 

2040  rem  the  insolation  calculations 

2050  rem  ************************************************* 

2060  rem  **  direct  normal  insolation  on  a  clear  day 
2070  idn  =  a  (mon)/(exp(b(mon)/sbeta)) 

2080  rem  **  total  horizontal  radiation  for  ktbar  day 
2090  ith  =  idn*(c(mon)+sbeta)*  pas(mon) 

2100  rem  **  sky  diffuse  for  ktbar  day 

2110  ids  =  i th*( 1 . 39-4. 03*ktb(mon)+5. 53*ktb(mon)~ 2-3. ll*ktb(mon)‘ 3) 
2120  rem  **  the  horizontal  direct  for  a  ktbar  day 
2130  idh  -  ith  -  ids 

2140  rem  **  direct  normal  for  a  ktbar  day 

2150  idn  =  idh  /  sbeta 

2160  rem  **  solar  direct  on  a  surface 

2170  ibt  =  idn*ctheta 

2180  rem  **  solar  diffuse  on  surface 

2190  idt  =  ids*((l+cti lt)/2)  +  i th*( ( 1-cti It))  *rho/2 

2200  rem  **  total  solar  on  surface 

2210  itt  =  ibt  +  idt 

2220  if  (itt  <=0)  then  itt  =  .00001 

2230  rem 

2240  rem  **********  calculate  new  ambient  temperature  ******* 

2250  kint  =  int(k) 

2260  tamb  =  (tavg(mon)+drange(mon)/2)-pcent(kin)*drange(mon) 

2270  rem 

2280  rem  **************************************************** 

2290  rem 

2300  rem  the  heat  transfer  calculations 

2310  rem  *************************************************** 

2320  qcond  =  ( tamb-ti n)/rtot 

2330  qconv  =  (tskin-tamb)*hconv 

2340  tsurf  =  (tamb+tskin)/2 

2350  qradout  =  eout*sigma*(tskin‘4  -  tsurf'4) 

2360  qradin  =  itt 
2370  gosub  5000 
2400  rem 

2410  rem  **  format  for  output  ** 

2420  time  =  kint  *  100 
2430  azm  =  azm/. 01745 
2440  psi  =  psi/. 01745 

2450  beta  =  atn( sbeta/sqr(-sbeta*sbeta+l))/. 01745 

2460  theta  =  (-atn(ctheta/sqr(-ctheta*ctheta+l))+c. 1416/2). 01745 

2470  rem 

2480  rem  **  output  ** 

2490  junk  =  k* 100  mon  100 

2500  if  junk  =  0  then  print  time,  tamb,  tskin 

2510  rem 

2520  rem  **  sum  the  radiation  for  the  day  ** 

2530  ibttot  =  ibttot  +  1bt*tstep 
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2540  idttot  =  idttot  +  idt*tstep 
2550  itttot  =  itttot  +  itt*tstep 
2560  tqradout  =  tqradout  +  qradout*step 
2570  tqcond  =  tqcond  +  qcond*step 
2580  tqconv  =  tqconv  +  qconv*step 

2590  rem 

2591  gosub  6000 

2592  rem 

2593  if  (tskin>tmax)  then  tmax  =  tskin 

2594  rem 
2600  next  j 

2610  rem  **  output  monthly  totals  ** 

2620  print  "qradin" ,  "qradout",  "qcond",  "qconv" 

2630  print  itttot,  tqradout,  tqcond,  tqconv 
2640  print 

2650  print 

2651  print  "The  maximum  skin  temperature  during  this  month  was", 
tmax 


2660  rem 
2670  next  mon 
2680  rem 

2685  next  direction 


2690  end 
2700  rem 

4000  rem******************************************************* 


4010 

rem- 

4020 

rem 

the  subroutine  to  calculate 

di rections 

4030 

rem 

4040 

rem 

orientation  m  e 

s 

w 

4050 

rem 

w  =  180  90 

0 

90 

4060 

rem 

w=  3.1416  1.5708 

0 

1.5708 

4070 

rem 

direction  =  0  for  north 

4080 

rem 

direction  =  1  for  east 

4090 

rem 

direction  =  2  for  west 

4100 

rem 

direction  =  3  for  south 

4110 

rem 

4120  if  (direction  =  0)  then  4130  else  4200 
4130  if  (k<12)  then  azm  =  (psi-3.1416) 

4140  if  (k>=12)  than  azm  =  (psi  +  3.1416) 

4150  rem 
4160  return 

4200  if  (direction  =  1)  than  4210  else  4300 
4210  azm  =  psi  -  w 
4220  return  4230  rem 

4300  if  (direction  =  2)  then  4310  else  4400 
4310  azm  =  psi  +  w 
4320  return 
4330  rem 

4400  if  (direction  =  3)  then  azm  =  psi 
4410  return 

5000  rem  ************************************************ 
5010  rem 

5020  rem  the  conduction  subroutine 
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************************************************ 


5030  rem 
5040  rem 
5050  condl  =  (tskin-tl)/rlayerl 
5060  storel  =  qradin-qradout-qconv-condl 
5070  cond2  =  (tl-t2)/rlayerl 
5080  store2  =  condl-cond2 
5090  cond3  =  ( t2-t3)/rl aterl 
5100  store3  =  cond2-cond3 
5110  cond4  =  (t3-t4)/rlayerl 
5120  store4  =  cond3-cond4 
5130  cond5  =  (t4-t5)/r1ayer2 
5140  store5  =  cond4-cond5 
5150  cond6  =  (t5-t6)/rlayer2 
5160  store6  =  cond5-cond6 
5170  cond7  =  (t6-t7)/rlayer2 
5180  store7  =  cond6-cond7 
5190  cond8  =  (t7-t8)/rlater2 
5200  store  =  cond7-cond8 
5210  cond9  =  (t8-t9)/rl 
5220  store9  =  cond8-cond9 
5230  condlO  =  (t9-tl0)/r2 
5240  storelO  =  cond9-condl0 
5250  return 

5260  rem  ************************************************ 
6000  rem 

6010  rem  the  temperature  calculations 

6020  rem  ************************************************* 

6030  dti  =  storel/(lll*shll*dl 1) 

6040  tskin  =  tskin  +dtl*tstep 
6050  dt2  =  store2/l 1 1 * sh 1 l*dll ) 

6060  tl  =  tl  +dt2*tstep 

6070  dt3  =  store3/( 111 *shl l*dll ) 

6080  t2  =  t2  +dt3*tstep 

6090  dt4  =  store4/( 1 1 l*shl l*dl 1 ) 

6100  t3  =  t3  +dt4*tstep 

6110  dt5  =  store5/(112*shl2*dl2) 

6120  t4  =  t4  +  dt4*tstep 

6130  dt6  =  store6/(112*shll2*dl2) 

6140  t5  =  t5  +  dt6*tstep 
6150  dt7  =  store7/( 112*shl2*dl2) 

6160  t6  =  t6  +  dt7*tstep 
6170  dt8  =  store8/( 112*shl2*dl2) 

6180  t7  =  t7  +  dt8*tstep 
6190  dt9  =  store9/( 113*shl3*dl3) 

6200  t8  =  t8  +  dt9*tstep 
6210  dtlO  =  storel0/( 114*shl4*) 

6220  t9  =  t9  +  dtl0*tstep 
6230  return 
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SAMPLE  OUTPUT 


MONTH  DAY  OF  THE  YEAR 

1  16 

********************* ************ 


THE  WALL  FACES  NORTH 


TIME 

TAMB 

TSKIN 

900.00 

28.28 

27.67 

1000.00 

31.01 

30.52 

1100.00 

34.10 

33.79 

1200.00 

37.01 

37.17 

1300.00 

39.20 

40.17 

1400.00 

40.65 

42.35 

1500.00 

41.20 

43.47 

1600.00 

40.65 

43.37 

THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS 

QRADIN  TOTAL  QRADOUT  TOTAL  QCOND  TOTAL 
166.3913  2. 415649E-03  129.9534 


MONTH  DAY  OF  THE  YEAR 

7  197 


******************************* 


THE  WALL  FACES  NORTH 


TIME 

TAMB 

TSKIN 

600.00 

64.60 

67.44 

700.00 

65.71 

72.41 

800.00 

67.72 

73.70 

900.00 

70.62 

76.23 

1000.00 

73.96 

79.98 

1100.00 

77.75 

84.04 

1200.00 

81.32 

88.15 

1300.00 

84.00 

91.36 

1400.00 

85.78 

94.70 

1500.00 

86.45 

96.43 

1600.00 

85.78 

96.85 

1700.00 

84.22 

98.20 

1800.00 

81.77 

99.56 

1900.00 

78.87 

95.49 

THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS 


QRADIN  TOTAL 
572.7928 


QRADOUT  TOTAL 
.228816 


QCOND  TOTAL 
265.8097 


43.61099 

QCONV  TOTAL 
36.43517 


99.55658 

QCONV  TOTAL 
306.7536 
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MONTH  DAY  OF  THE  YEAR 

1  16 


********************************** 


THE  WALL  FACES  SOUTH 


TIME 

TAMB 

TSKIN 

900.00 

28.28 

38.22 

1000.00 

31.01 

48.74 

1100.00 

34.10 

58.51 

1200.00 

37.01 

67.05 

1300.00 

39.20 

73.55 

1400.00 

40.65 

77.27 

1500.00 

41.20 

77.69 

1600.00 

40.65 

73.91 

THE  MAXIMUM 

SKIN  TEMPERATURE  Dl 

QRADIN  TOTAL 

QRADOUT  TOTAL 

929.9351 

.1486741 

MONTH 

DAY 

OF  THE  Yl 

7 

197 

******************************: 

THE  WALL  FACES  SOUTH 

TIME 

TAMB 

TSKIN 

600.00 

64.60 

65.01 

700.00 

65.71 

67.23 

800.00 

67.72 

70.22 

900.00 

70.62 

76.23 

1000.00 

73.96 

84.24 

1100.00 

77.75 

92.48 

1200.00 

81.32 

99.99 

1300.00 

84.00 

105.78 

1400.00 

85.78 

108.99 

1500.00 

86.45 

109.19 

1600.00 

85.78 

106.24 

1700.00 

84.22 

102.76 

1800.00 

81.77 

98.88 

1900.00 

78.87 

94.33 

qCOND  TOTAL 
362.7757 


THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS 


QRADIN  TOTAL 
736.0322 


QRADOUT  TOTAL 
.4497645 


QCOND  TOTAL 
271.1892 


78.05954 

QCONV  TOTAL 
567.0099 


109.6272 

QCONV  TOTAL 
464.3918 
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MONTH  DAY  OF  THE  YEAR 

1  16 


********************************** 


THE  WALL  FACES  WEST 


TIME 

TAME 

TSKIN 

900.00 

28.28 

27.67 

1000.00 

31.01 

30.52 

1100.00 

34.10 

33.79 

1200.00 

37.01 

37.17 

1300.00 

39.20 

43.48 

1400.00 

40.65 

51.28 

1500.00 

41.20 

58.26 

1600.00 

40.65 

61.70 

THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS 

QRADIN  TOTAL  QRADOUT  TOTAL  QCOND  TOTAL 
402.7702  2.352482E-02  248.3496 


MONTH  DAY  OF  THE  YEAR 

7  197 


********************************* 


* 


THE  WALL  FACES  WEST 


TIME 

TAME 

TSKIN 

600.00 

64.60 

65.01 

700.00 

65.71 

67.23 

800.00 

67.72 

70.22 

900.00 

70.62 

73.93 

1000.00 

73.96 

78. 19 

1100.00 

77.75 

82.63 

1200.00 

81.32 

87.04 

1300.00 

84.00 

95.78 

1400.00 

85.78 

107.13 

1500.00 

86.45 

118.50 

1600.00 

85.78 

128.23 

1700.00 

84.22 

134.57 

1800.00 

81.77 

135.11 

1900.00 

78.87 

121.55 

THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS 


QRADIN  TOTAL 
1149.943 


QRADOUT  TOTAL 
1.006868 


QCOND  TOTAL 
504.2972 


61.70818 

QCONV  TOTAL 
154.3973 


135.6975 

QCONV  TOTAL 
644.6388 
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MONTH  DAY  OF  THE  YEAR 

1  16 


****** **************************** 


THE  WALL  FACES  EAST 


TIME 

TAMB 

TSKIN 

900.00 

28.28 

38.97 

1000.00 

31.01 

44.87 

1100.00 

34.10 

47.21 

1200.00 

37.01 

46.74 

1300.00 

39.20 

47.08 

1400.00 

40.65 

47.76 

1500.00 

41.20 

47.73 

1600.00 

40.65 

46.73 

THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS  47.98488 

QRADIN  TOTAL  QRADOUT  TOTAL  QCOND  TOTAL  QCONV  TOTAL 
364.2411  1.756502E-02  161 . 41  IE  202.8114 


MONTH  DAY  OF  THE  YEAR 

7  197 


********************************** 


THE  WALL  FACES  EAST 


TIME 

600.00 

700.00 

800.00 

900.00 

1000.00 

1100.00 

1200.00 

1300.00 

1400.00 

1500.00 

1600.00 

1700.00 

1800.00 

1900.00 


TAMB 

64.60 

65.71 

67.72 
70.62 
73.96 
77.75 
81.32 
84.00 
85.78 
86.45 
85.78 
84.22 
81.77 
78.87 


TSKIN 

72.53 

91.53 

102.99 
110.27 
113.73 
113.44 
109.80 
108.00 
107.37 
106.41 
104.71 

101.99 
98.31 
93.89 


THE  MAXIMUM  SKIN  TEMPERATURE  DURING  THIS  MONTH  WAS  114.4237 

QRADIN  TOTAL  QRADOUT  TOTAL  QCOND  TOTAL  QCONV  TOTAL 
1100.46  .876272  274.3008  825.2806 
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Appendix  C 


Dimensions  and  Locations  of  Buildings 
Included  in  Heat  Rejection  Analyses 


All  buildings  along  that  portion  of  Wisconsin  Avenue  passing  to  the 
southwest  of  the  Naval  Observatory  were  included  in  these  analyses.  Each 
is  identified  by  a  number  on  Figure  C-l.  Sizes,  HVAC  loads,  and  loca¬ 
tions  of  the  buildings  are  listed  on  Table  C-l.  Recall  that  it  is  being 
assumed  that  all  buildings  are  square.  Ten  feet  (three  meters)  of  height 
is  used  for  each  floor.  The  numbers  in  the  "Location"  column  refer  to 
the  sectioning  of  Wisconsin  Avenue  to  simplify  the  analyses.  This  was 
discussed  previously;  see  Figure  22. 

The  "typical  building"  at  the  bottom  of  Table  C-l  is  the  baseline 
building  referred  to  in  Table  5  and  Figures  24  and  25.  It  is  assumed  to 
be  located  directly  west  of  the  26-inch  refracting  telescope,  and  desig¬ 
nated  as  No.  19  on  Figure  C-l. 


C-l 


Figure  C-1 .  Schematic  of  the  Naval  Observatory  and  surroundings  showing  the 
buildings  included  in  the  heat  rejection  calculations. 
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33  Site  data  and  systems  integration  (energy  resource  data. 

energy  consumption  data.  Integrating  energy  systems) 

34  ENVIRONMENTAL  PROTECTION 

35  Hazardous  waste  minimization 

36  Restoration  of  Installations  (hazardous  waste) 

37  Waste  water  management  and  sanitary  engineering 

38  Oil  pollution  removal  and  recovery 

39  Air  pollution 

44  OCEAN  ENGINEERING 

45  Seafloor  soils  and  foundations 

46  Seafloor  construction  systems  and  operations  (including 

diver  and  manipulator  tools) 

47  Undersea  structures  and  materials 

48  Anchors  and  moorings 
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50  Pressure  vessel  facilities 

51  Physical  environment  (including  site  surveying) 

52  Ocean-based  concrete  structures 
54  Undersea  cable  dynamics 
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